nTlc  ntE  COPY 


Best  Available  Copy 

AD- A 167  363 

Report  021779-2-T 


SIGNAL  DETECTION  AND  PARAMETER 
IDENTIFICATION  WITH  APPLICATION 
TO  UNDERWATER  ACOUSTICS 


Fuad  Imran  Khan 


COMMUNICATIONS  AND  SIGNAL  PROCESSING  LABORATORY 

Department  of  Electrical  Engineering  and  Computer  Science 
The  University  of  Michigan 


Ann  Arbor,  Michigan  48109 


February  1986 


Technical  Report  No.  233 
Approved  for  public  release;  distribution  unlimited 


S 


DTIC 

ELECTE 
MAY  0  7  1986 


Prepared  for 

OFFICE  OF  NAVAL  RESEARCH 

Department  of  the  Navy 
Arlington,  Virginia  22217 

86  5  5  053 

Best  Available  Copy 


.N 

S 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 

la.  REPORT  SECURITY  CLASSIFICATION 

1b.  RESTRICTIVE  MARKINGS 

Unclassified 

None 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  DISTRIBUTION /AVAILABtUTY  OF  REPORT 

Approved  for  puhlie  relAasof  di  st-ri  hn-hion 

2b.  DECLASSIFICATION /downgrading  SCHEDULE 

unlimited. 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

021779-2-T 

S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

1 

j 

6c  ADDRESS  (Cty.  Staf,  amiZI^Codt) 
The  University  of  Michigan 

Ann  Arbor,  Michigan  48109 


6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  apfilktbh)  Office  of  Naval  Research 

Code  425UA 


7b.  ADDRESS  (Oty,  State,  and  Z/A  Coda) 
800  North  Quincy  Street 

Arlington,  Virginia  22217 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZAHON 


Be  ADDRESS  (Oty,  State,  and  2/P  Code) 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(tfappP-aM) 

Contract  No.,  N00014-84-K-0066 


to.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT  TASK 

ELEMENT  NO.  NO.  NO. 

NR  386-941 


1 1 .  TITLE  (Includa  Saturit^  Classitkation) 

SIGNAL  DETECTION  AND  PARAMETER  IDENTIFICATION 

WITH  APPLICATION  TO  UNDERWATER  ACOUSTICS 

12.  PEMONAL  AUTHOR(S) 

Ktian,  Fuad  Imran 

13a.  TYPE  pF  REPORT 

Technical  Report 

13b.  TIME^CpVEREO  „ 

FROM  TO  2/a6 

14.  DATE  OF  REPORT  (Y*ar,  Month,  Day) 
February  1986 

IS.  PA^E  COUNT 

COSATt  CODES 


GROUP  i  SUB-GROUP 


18.  SUBJECT  TERMS  (Cont/nut  on  rtMisa  it  imaiisty  and  kHirtify  bjf  bhtk  nmbat) 
Passive  Detection  Multipath  Acoustic  Propagation 

Underwater  Acoustics 


19.  ABST^CT  (Cont/ntff  on  revenw  it  nectsaty  and  kkndfy  by  btodt  nwnbati 

A  problem  of  passive  detection  of  underwater  moving  acoustic  sources  is  formulated. 

It  is  found  that  solving  the  exact  problem,  when  a  multipath  acoustic  propagation  environ¬ 
ment  is  assumed,  is  fairly  difficult.  Keeping  the  multipath  assumption  intact,  we  develop 
two  abstract  problems  that  are  solvable  and  yet  are  sufficiently  related  to  reality  to 
allow  us  to  form  guidelines  and  rules  of  thumb  for  the  passive  detection  problem.  Multipatli 
propagation  from  a  moving  acoustic  source  is  characterized  by  doppler  and  delay  parameters 
on  each  path.  We  assume  that  sound  is  received  at  two  remote  receivers.  If  there  are  M 
patiis  to  one  receiver  and  N  paths  to  the  other  receiver,  the  normalized  crosscorrelation 
of  the  receptions  will  exhibit  M.N  peaks. 

We  use  tlie  concept  of  Extended  M-orthogonal  signals  to  unify  t)ie  information  on  the 
ambiguity  plane  and  to  find  how  many  of  tiie  peaks  on  the  ambiguity  plane  need  to  be 
localized  for  acceptable  detector  perform€mce.  -pi-  - 

The  M.N  peaks  are  generated  by  M+N  independent  Mts  of  doppler  and  delay  parametei^s^^^ 


20.  D'STR'BUTION/AVAILAIILITY  OF  abstract  21.  ABST/UCT  SECURITY  CLASSIFICATION  ] 

□  UNCLASSIFIEOAJNLIMITED  ED  SAME  AS  RPT.  □  OTIC  USERS  pnclassified 


22*.  NAME  OF  RESPONSIBLE  INDIVIDUAL  |22b.TElEPHONE(lnc/!K)»  Ana 

Carol  S.  Van  Aken  _ _  I  [313)  764-5220 


OOFORM ia3.B«MAR  B3 APR t<«t«>nm.ybtu«<} until tKhidiUd.  SECURITY  ur  .m. 

All  other  editiont  art  obsol«U.  - 

UNCLASSIFIED 

 i 


22c.  OFFICE  SYMBOL 


UNCLASSIFIED 

HCMWTY  CkAWnCATlOW  Of  TWt  >»Aai 

19.  Abstract  (cent) . 

Periodic  Rcindom  Sequences  are  used  to  address  the  question: 

Should  the  detector  treat  the  M.N  peeiks  as  independent  entities  or  as  being 
generated  from  M+N  independent  parameter  sets?  i.e.  Should  the  detector  have  M.N 
degrees  of  freedom  or  M+N  degrees  of  freedom? 


ACKNOWLEDGEMENTS 


I  would  like  to  thank  the  members  of  my  dissertation 
committee,  Profs.  Neuhoff,  Perakis,  Root  and  Stark,  Their 
suggestions  and  guidance  made  this  a  better  thesis.  In 
particular  I  would  like  to  thank  Prof.  Birdsall,  who  over 
the  years  has  been  both  a  teacher  and  a  friend.  Without  his 
guidance,  insight  and  inspiration  my  graduate  studies  would 
have  been  incomplete. 

I  would  like  to  thank  the  staff  of  Cooley  Electronics 
Lab.  for  their  support,  in  particular  I  would  like  to 
mention  Dr.  Kurt  Metzger  for  several  useful  discussions  and 
for  help  with  the  word  processor,  Ms.  Linda  Russell  for  her 
help  with  FORM  and  the  word  processor.  Mr.  John  Campbell 
for  his  help  with  programming,  Ms,  Carol  Van  Aken  for  her 
administrative  assistance.  Mr.  Wayne  V,  Wald  for  his 
assistance  with  the  figures.  I  am  grateful  to  my  colleagues 
and  friends  Matt  and  Deb  for  their  encouragement  and  help. 

My  friends  Babar,  J.B.,  Naveed«  Alan  and  Firdous 
encouraged  and  inspired  me  to  excel,  I  would  like  to  thank 
them  for  this  and  much  more.  I  would  like  to  thank  Feisal 
for  his  help  with  the  figures.  I  would  like  to  thank  Aliya 
for  her  encouragement  and  humor.  Last  but  not  least  1  wish 
to  thank  my  dear  friend  Farzana  for  her  encouragement, 
support  and  understanding  in  both  good  tioies  and  bad. 

Finally  I  would  like  to  thank  my  parents  and  my  sister 
for  providing  me  with  the  opportunity  of  studying  in  the 
United  States  and  for  their  constant  support,  encouragement 
and  love. 


The  support  of  the  Office  of 
42SUA}  is  gratefully  acknowledged. 


•  *  * 
111 


Naval  Research  (Code 


TABLE  OP  CONTENTS 


Page 

ACKNOWLEDGEMENTS  i i i 

LIST  OF  TABLES  vii 

LIST  OF  FIGURES  ix 

LIST  OF  APPENDICES  xi 

LIST  OF  SYMBOLS  AND  ABBREVIATIONS  xiii 

CHAPTER 

I .  INTRODUCTION  1 

II.  REVIEW  OF  DETECTION  THEORY  7 

2.1  The  Basic  Detection  Problem  7 

2.2  Optimum  Receivers,  Detection  Criteria 

and  Decision  Rules  9 

2.3  Simple  and  Composite  Signal  Hypothesis  11 

2.4  Performance  Evaluation  13 

2.5  Two  Examples  14 

2.5.1  Known  Scalar  Signal  in  Gaussian  Noise  14 

2.5.2  Gaussian  Signal  in  Gaussian  Noise  16 

III.  PROBLEM  FORMULATION  21 

3.1  The  Underwater  Acoustic  Channel  21 

3.2  Receiver  Geometry  and  Hypothesis 

Formulation  23 

3.3  Post  Reception  Processing  for  Detection  25 

IV.  EXTENDED  M-ORTHOGONAL  SIGNALS  AND  THE 

FIRST  ABSTRACTION  33 

4.1  Abstraction  Formulation  33 

4.2  The  Orthogonal'-Signal-'Based  Receiver  36 

4.2.1  The  M-Orthogonal  Signals  Evaluation  37 

4.2.2  The  Extended  M-Orthogonal  Signals 

Evaluation  37 

4.3  Receiver  Performance  Evaluation  39 

4.3.1  Approximation  to  the  Detectability 

Index  42 

4.3.2  Formal  Evaluation  of  the  pdf  of  £(y|K}  44 

4.3.3  Discrete  Fourier  Transform  Approach  47 

4.3.4  Numerical  Integration  Methods  for 

Evaluating  the  R.O.C.  48 

4.3.5  Monte-Carlo  Methods  for  Evaluating 

the  R.O.C.  54 

4.4  Summary  57 


V 


V.  PERIODIC  RANDOM  SEQUENCES  AND  THE  SECOND 

ABSTRACTION  59 

5.1  Periodic  Random  Sequences  59 

5.2  Detection  of  one  PRS  of  Known  Period  in 

Gaussian  Noise  60 

5.2.1  Properties  of  Circulant  Matrices  62 

5.2.2  Eigenvalue  Eigenvector  Decomposotion 

of  R  64 

5.2.3  The  Optimum  Detector  65 

5.2.4  Performance  Evaluation  of  the 

Optimum  Detector  72 

5.2.5  Detection  of  PRS*s  with  Unknown 

Period  84 

5.3  Detection  of  two  or  more  PRS's  in  Noise  97 

5.3.1  Detection  of  2  or  more  PRS's  of 

Known  Periods  97 

5.3.2  Detection  of  k  of  M  PRS's  in  Noise  114 

5.3.3  Detection  of  Arbitrary  PRS  in  non 

White  Gaussian  Noise  121 

5.4  The  Abstract  Problem  Revisited  124 

VI.  SUMMARY  AND  CONCLUSIONS  127 

APPENDICES  131 


REFERENCES 


145 


LIST  OF  TABLES 


Table 

4.1 


Title 


Table  of  normal  detectabiltiy  for 
several  K  values,  for  signal  power  S"!. 

Statistics  of  the  random  variable  W|^ 
for  K-0, 1,2, 4, 8, 16, 32, 64. 


4.2 


LIST  OP  FIGURES 


Fiqure 

Title 

Page 

2.1 

Block  diagram  of  the  basic  detection 
problem 

8 

2.2 

Normal  R.O.C*  curves. 

17 

3.1 

Receiver  Geometry. 

24 

3.2 

Block  diagram  of  the  post  reception 
processing. 

27 

4.1 

R.O.C.  curves  for  the  standard 

Extended  M-Orthogonal  signals. 

55 

4.2 

R.O.C.  curves  for  the  modified 

Extended  M'-Orthogonal  signals. 

56 

5.1 

Block  diagram  of  the  Frequency  Selective 
Energy  Detector. 

68 

5.2a 

Block  diagram  of  the  Circulating  Average 
Energy  Detector  Signal  Processor  vith 
Parameters  P^,  Li  and  Ai. 

71 

5.2b 

Block  diagram  of  the  Circulating  Average 
Energy  Detector. 

73 

5.3 

Block  diagram  of  the  Energy  Detector. 

77 

5.4a 

Averaged  R.O.C.  curves  for  the  CAED. 

79 

5.4b 

R.O.C.  curves  for  the  ED. 

80 

5.5 

Block  diagram  of  the  optimum  detector 
when  the  signal  period  is  equally 
likely  to  be  one  of  M  periods. 

87 

5.6 

Averaged  R.O.C.  curves  for  the  optimum 
detector  vith  M<*10. 

88 

5.7 

Block  diagram  of  the  estimator  detector 
when  the  signal  period  is  equally 
likely  to  be  one  of  M  periods. 

89 

5.6 

Averaged  R.O.C.  curves  for  the  estimator 
detector  vith  N-10* 

90 

5.9 

R.O.C.  curves  for  the  energy  detector. 

91 

5,10 

Average  Prob.  of  Error  in  estimating 

ix 


pucvious  PAoe 

IS  BLANK 


9 


93 


signal  period  vs.  S/N  ratio. 

5.  It  Block  diagram  of  the  estimator  estimator.  98 

5.12  R.O.C.  curves  for  the  optimum  detector 

for  2  PRS's  of  known  periods  in  noise.  105 

5.13  R.O.C.  curves  for  the  sum  detector  for 

2  PRS's  of  known  periods  in  noise.  106 

5.14  R.O.C.  curves  for  the  energy  detector.  107 

5.15  Averaged  R.O.C.  curves  for  the  quasi 

optimum  detector  for  N>2.  110 

5.16  Averaged  R.O.C.  curves  for  the  quasi 

sum  detector  for  N>2.  Ill 

5.17  R.O.C.  curves  for  the  energy  detector.  112 

5.18  Averaged  R.O.C.  curves  for  the  quasi 

optimum  detector  for  M«>10  and  ke2.  118 

5.19  Averaged  R.O.C.  curves  for  the  estimator 

detector  for  M«10  and  k«2.  119 

5.20  Averaged  R.O.C.  curves  for  the  max 

detector  for  M«10  and  k«2.  120 


X 


LIST  OF  APPENDICES 


Appendix 

A.  Monte  Carlo  Simulation  Methods. 

B.  Derivation  of  the  Normal  Detectability 
of  the  PSED. 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


BNC 

CAEO 


DPT 

Do 

Di 

d 

d’ 

dB 

ED 

CMOS 

B(.) 

FSEO 

£,,(.) 


Ho 

Hi 


i.i.d. 


L 


Power  level  of  Periodic  Random  Sequence  with 
period  P. 

Bivariate  Normalized  Crosscorrelation 
Circulating  Average  Energy  Detector 
Generic  nxn  circulant  matrix 

Normalized  complex  attenuation  coefficient 
of  the  nth  path 

Discrete  Fourier  Transform 

Decision  in  favor  of  the  noise  alone  hypothesis 

Decision  in  favor  of  the  signal'<*noise  hypothesis 

Normal  detectability 

Square  root  of  the  normal  detectability 

Decibel 

Energy  Detector 

Extended  M-^Orthogonal  signals 

Expected  value 

Frequency  Selective  Energy  Detector 

Probability  d  nsity  function  of  random 
variable  x 

Conditional  Probability  density  function 
Noise  alone  hypothesis 

Signal'^'Noise  hyi>othesis 

LxL  identity  matrix 

independent  and  identically  distributed 
Number  of  samples  in  observation 

pnevious  paoc 

IS  RLANK 

xiii 


A(.) 

£(.l.) 

MOS 

N(v,o  ) 
n(.) 


n 


*^FA 

PRS 


p  •  d  •  £ « 
Pi(.) 


pj(.y 

R 

R(.) 


R.O.C. 


r.v» 


s( . ) 


5 


U(0,1) 

Y(.) 


Number  of  integer  periods  of  PRS  with  period 
P^  in  the  observation 

Least  Common  Multiple 
Likelihood  ratio 
Conditional  likelihood  ratio 
M-Orthogonal  signals 

2 

Normal  density  with  mean  u  and  variance  a 
Noise 

Noise  vector 

Probability  of  Detection 

Probability  of  False  Alarm 

Periodic  Random  Sequence 
probability  density  function 
PRS  with  period  Pj 

One  period  of  PRS  p^C.) 

Autocorrelation  matrix 
Autocorrelation  function 
Receiver  Operating  Characteristic 
random  variable 
Signal 

Signal  vector 
Estimate  of  s 

Unitary  matrix  of  the  eigenvectors  of  a 
nxn  circulant  matrix 

Uniform  density  over  [0,1) 

DFT  of  the  sequence  y(.} 


xiv 


y  Vector  of  the  sequence  Y(.) 

y(.)  Observation 

^  Observation  vector 

z(.)  Log-likelihood  ratio 

2(,|.)  Conditional  log-likelihood  ratio 

z(.)  Quasi  log-likelihood  ratio 

g  Differential  doppler  correction  parameter 

$(.)  Delta  function 

$(.)  Normal  distribution  function 

(Ji  ( . )  Characteristic  function  of  random  variable  x 

A 

r(.)  Gamma  function 

r(.,.)  Incomplete  gamma  function 

y(.,.)  Incomplete  gamma  function 

Bivariate  Normalized  Crosscorrelation  function 
with  an  observation  interval  of  T  seconds 

y(.,.)  Bivariate  Normalized  Crosscorrelation  function 

with  an  infinite  observation  interval 


Y  Observation  pattern  in  the  discrete  ambiguity 

plane 

X  Vector  associated  with  the  observation 

pattern  in  the  discrete  ambiguity  plane 

2 

0^  Variance  of  random  variable  x 

T  Differential  time  delay  correction  parameter 

Tyj(t)  Time  delay  on  tie  nth  path  at  time  t 

Time  delay  on  the  nth  path  at  time  0 

Time  delay  derivative  of  the  nth  path 

Diagonal  matrix  of  the  eigenvalues  of  an 


nxn  circulant  matrix 


Ambiguity  function  with  an  observation 
interval  of  T  seconds 

Ambiguity  function  with  an  infinite 
observation  interval 


CHAPTER  I 


INTRODUCTION 

Passive  Sonar  Detection  and  the  associated  problems  of 
doppler  and  time  delay  estimation  are  problems  of  continuing 
research  interest  in  the  underwater  acoustics  community 
[1, 2,3,4]*  While  passive  systems  have  the  disadvantage  that 
they  can  not  control  signal  energy,  their  anonymity  and 
overall  cost  dictate  their  feasibility  especially  in 
surveillance  systems  [2], 

A  passive  sonar  system  utilizes  an  array  of  two  or  more 
receivers  (hydrophones).  Two  receivers  are  sufficient  to 
estimate  the  bearing  angle  of  the  source  in  the  plane  of  the 
receivers  [2],  However  if  more  information  is  required 
about  the  source,  such  as  localization  in  three  dimensions 
and  estimation  of  source  velocity,  more  receivers  are 
required.  The  accuracy  of  the  estimates  is  a  function  of 
the  receiver  geometry  in  three  dimensions  as  well  as 
measurement  accuracy.  Throughout  this  thesis  we  will  assume 
that  we  have  two  receivers.  This  is  a  common  assumption 
made  by  several  authors  [1,2,5].  It  is  not  too  difficult  to 
extend  the  analysis  to  more  than  two  receivers  if  the 
receivers  are  considered  in  pairs* 

We  assume  a  multipath  channel  for  sound  propagation  to 
each  of  the  receivers.  This  is  in  contrast  to  the 
assumption  usually  made;  i.e.  there  is  only  one  path  to 
each  receiver  [2,6].  Associated  with  each  path  in  the 
channel  is  a  complex  attenuation  (or  gain)  coefficient,  a 
doppler  parameter  due  to  source  motion  and  time  delay  in 
propagation  from  the  source  to  the  receiver.  The  details  of 
the  multipath  channel  are  given  in  Chapter  ill.  The 
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implications  of  the  multipath  channel  assumption  will  become 
clear  as  we  proceed. 

We  will  use  the  Bivariate  Normalized  Crosscorrelation 
(BNC)  function  as  the  observable  for  our  detector. 
Bivariate  refers  to  the  presence  of  two  search  variables. 
We  define  the  BNC  function,  denoted  by  as  follows 


rT(s-T) 

(?  /t  |yi(t)|*dt)’/^(i |y2(t)|'dt)’/2 

where  yi(t)  and  y2(t)  denote  the  receptions  at  the  two 
rec**lve-3,  /.j,  signifies  an  integral  over  an  interval  of 
length  T.  We  will  be  interested  in  the  intervals  tO,T)  and 
[-T/2,T/2).  In  general  is  complex.  A  simple 

replication  of  the  Schwarz  inequality  shows  that 
IYt<B^t)|  <  1.  We  will  also  need  to  consider  the  ambiguity 
function  i|),p(B,T)  defined  as 


♦,,(8,t)  |YT(8.t)|'  <1-2) 

The  properties  of  the  ^NC  function  and  the  ambigutiy 
function  are  extensively  documented  in  the  literature  [7]. 
We  will  also  have  occasion  to  consider  tha  BNC  and  the 
ambiguity  functions  in  the  limit  T  We  define 

Y(B»t)  ■  lim  Y«rU»T)  (1.3a) 

4»(B,t)  »  lim  (i.3b) 

T^cd  * 


The  search  variables  0,  t  seek  to  O’atch  differential 
doppler  and  differential  delay  between  paths  to  the  two 
receivers.  Two  input  cot  relators,  univariate  and 
multvariate,  have  been  used  extensively  in  the  detection 


3 


theory  literature  to  form  the  observation  statistic  [1,8,9]. 

A  detector  that  is  based  on  the  assumption  that  there 
is  a  single  path  to  each  receiver  seeks  to  identify  a  single 
maximum  in  the  ambiguity  diagram.  The  maximum  is  related  to 
the  mean  differential  delay  and  the  mean  differential 
doppler.  For  detection  purposes  the  height  of  the  maxmimum 
is  compared  to  a  pre-selected  threshold.  Detection  is 
indicated  if  the  height  of  the  maximum  is  greater  than  the 
threshold.  However  a  detector  that  is  based  on  a  multipath 
channel  model  seeks  to  identify  all  the  resolvable  path 
pairs.  So,  if  there  are  M  paths  to  one  receiver  and  N  paths 
to  the  other  receiver  and  all  the  path  pairs  are  resolvable, 
the  detector  based  on  the  multipath  channel  model  will 
identify  MN  peaks  (local  maxima)  in  the  ambiguity  plane.  To 
resolve  the  peaks  (path  pairs),  we  need  integration  times 
that  are  long  [10,11,12].  More  is  said  about  integration 
times  in  Chapter  III. 

From  the  above  it  follows  that  if  the  sound  propagation 
is  by  a  multipath  channel,  then  a  detector  based  on  the 
multipath  channel  model  has  more  apriori  information 
available  to  it  as  compared  to  a  detector  based  on  the 
single  path  channel  model.  When  sound  propagation  is  by  a 
multipath  channel,  we  expect  the  performance  of  the 
"multipath  channel  detector"  to  be  superior  to  the 
performance  of  the  "single  path  channel  detector."  In 
making  a  decision  an  optimum  detector  utilizes  all  the 
information  available  about  the  signal  (or  source)  in  the 
observation.  The  primary  goal  of  this  thesis  is  not  to  find 
exact  analytical  solutions  but  to  establish  guidelines  and 
rules  of  thumb  about  how  best  to  unify  the  information 
available  in  the  ambiguity  plane  as  a  detection  statistic. 

Outline  of  the  Thesis 

Chapter  II  is  a  brief  review  of  fixed  time  detection 
theory,  it  is  included  for  completeness  and  for  establishing 


notation  and  convention. 

Chapter  III  starts  off  with  a  description  of  the 
multipath  channel  model.  After  that  the  problem  is 
formulated  and  developed  analytically.  It  is  found  that 
analytical  solutions  are  non-trivial  and  perhaps  not 
possible.  This  lays  the  foundation  for  the  abstract 
problems  considered  in  Chapters  IV  and  V. 

The  object  of  formulating  abstract  problems  is  to  study 
problems  that  aro  sufficiently  related  to  reality  and  yet 
are  solvable.  The  information  gained  from  the  study  of  the 
abstract  problems  is  then  utilized  to  form  guidelines  and 
rules  of  thumb  for  the  exact  problem. 

In  Chapter  IV  we  study  the  first  of  two  abstractions. 
Here  we  assume  a  discrete  ambiguity  plane  and  use  the 
"Extended  M-Orthogonal  Signals"  approach  to 

1)  Unify  the  information  in  the  ambiguity  plane 

2)  Find  how  many  peaks  (path  pairs)  need  to  be  matched 
for  good  or  acceptable  performance. 

The  MN  peaks  in  the-  ambiguity  diagram  are  generated  by 
M-«-N  independent  sets  of  doppler  and  delay  parameters. 
Chapter  V  uses  the  notion  of  "Periodic  Random  Sequences”  to 
address  the  following  questions 

1)  Should  the  detector  treat  the  MN  peaks  as 
independent  entities  or  as  being  generated  from  M't-N 
independent  parameter  sets,  i.e.  should  the  detector  have 
MN  degrees  of  freedom  or  M-^-N  degrees  of  freedom? 

2)  Does  one  approach  offer  a  significant  improvement  in 
performance  over  the  other? 

In  Chapter  VI  we  summarize  and  highlight  the  findings 
of  the  study.  Some  suggestions  for  future  work  are  also 
made. 

Some  liberty  has  been  taken  with  the  notations,  however 
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the  notations  are  consistent  and  the  author  hopes  they  are 
clear  from  the  context. 


PREVIOUS  PAGE 
IS  BLANK 


CHAPTER  II 

REVIEW  OF  DETECTION  THEORY 

The  purpose  of  this  chapter  is  to  review  some  of  the 
results  of  detection  theory  that  will  be  used  in  later 
chapters.  This  will  also  serve  to  establish  notation  and 
convention.  Since  most  of  the  results  presented  here  are 
well  known  [13,14,15,16,17]  detailed  derivations  are  not 
emphasized.  For  the  most  part,  we  will  follow  a  Bayesian 
approach  throughout  the  thesis.  By  this  we  mean  that  signal 
and  noise  parameters  that  are  not  known  exactly  can  be 
modelled  as  random  variables  with  known  probability  density 
functions  (p.d.f.'s).  The  p.d.f.'s  are  chosen  to  reflect 
the  observer's  state  of  objective  and  subjective  knowledge 
about  the  unknown  parameters. 

2.1  The  Basic  Detection  Problem 

A  block  diagram  of  the  basic  problem  in  fixed  time 
detection  theory  is  given  in  fig.  2.1.  The  observation  y(t) 
belongs  either  to  the  noise  alone  hypothesis  denoted  Ho  or 
to  the  signal-t-noise  hypothesis  denoted  Hi. 

Ho  :  y(t)  ■  n(t)  (2.1a) 

Hi  :  y(t)  ■  n(t)  +  8(t)  (2.1b) 

0  <  t  <  T 

Where  s(t)  is  a  known  signal  or  a  random  process.  It  is 

assumed  that  the  hypothesis  in  effect  does  not  change  over 


Block  diagram  of  the  basic  detection  problem 
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the  observation  interval.  In  other  words  the  hypotheses  H© 
and  Hi  are  mutually  exclusive. 

Often  for  analytical  convenience  random  processes  are 
expressed  as  random  vectors.  The  elements  of  the  random 
vector  may  be  obtained  using  the  Shannon  sampling  theorem 
[16,  pp.  29-301  or  by  using  a  truncated  form  of  some  series 
expansion  such  as  the  Karhunen-Loeve  expansion  [16,  pp.  54- 
74].  With  vector  notation  the  observation  is  written  as 

Ho  :  Z  ■  n  (2.2a) 

Hi  :  y  ■  n  +  s  (2.2b) 

where  y,  n  and  s  are  column  vectors  with  a  fixed  number  of 
elements,  say  L.  The  observation  under  each  hypothesis  is 
then  characterized  by  the  joint  p.d.f.'s  fo(y)  and  fi(y) 
where 

^  ^ (y 0 #y i * •  •  * ^ I )  t  i*o,i  (2.3) 

and  if  given  H>  the  elements  of  y  are  independent  and 
identically  distributed  (i.i.d.)  we  have 


L-1 

f4(y)  •  n  f(yu|Hi)  (2.4) 

k«0 

2.2  Optimum  Receivers  Detection  Criteria  and  Decision  Rules 

Based  on  the  observation  and  a  detection  criterion  the 
processor  computes  a  detection  statistic  a(y).  The  decision 
device  makes  a  binary  decision  Dq  or  Di,  according  to  a 
decision  rule  g(y)  and  the  value  of  a(y).  Do  implies  a 
decision  by  the  receiver  that  y  is  from  Ho;  Di  implies  a 
decision  by  the  receiver  that  y  is  from  Hi.  There  are  four 
possible  hypothesis-decision  pairs 
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1)  Decide  Do  when  Hq  is  true  :  Correct  Rejection  (CR) 

2)  Decide  Di  when  Hq  is  true  :  False  Alarm  (FA) 

3)  Decide  Dq  when  Hj  is  true  :  Miss  (M)  (2.5) 

4)  Decide  Di  when  Hi  is  true  :  Detection  (D) 

The  conditional  probabilities  associated  with  the 
hypothesis-decision  pairs  are 


1)  Probability  of  Correct  Rejection  ^  «  P(Do|Ho) 

2)  Probability  of  False  Alarm  ^  P^^  >  P(Di|Ho) 

3)  Probability  of  Miss  ^  P,4  -  P(Do|Hi)  (2.6) 

4)  Probability  of  Detection  ^  Pq  ■  P(Di|Hi) 

The  design  of  the  optimum  receiver  is  based  on  the 
conditional  probabilites  associated  with  the  hypothesis- 
decision  pairs. 

Commonly  used  detection  criteria  are  the  Bayes 
criterion,  the  Minimax  criterion  and  the  Neyman-Pearson 
criterion.  Based  on  available  knowledge  all  criteria  seek 
to  minimize  the  cost  associated  with  making  a  decision.  It 
has  been  shown  by  Birdsall  [16]  that  for  the  above  criteria 
and  any  other  criteria  that  consider  correct  decisions 
"good"  and  incorrect  decisions  "bad"  an  optimum  detection 
statistic  is  the  likelihood  ratio  £(^) 


A(2>  ^  <2.7) 

Upon  making  an  observation  the  optimum  receiver 
computes  the  likelihood  ratio  £(y;)»  The  decision  device 
then  compares  £(^)  with  a  pre-assigned  threshold  c  (0<c<<») 
and  decides  Di  if  £(^)  >  c  and  Do  if  £(y)  <  c.  If  £(^)  ■>  c 
then  with  probability  6  decision  Di  is  made,  this  is  called 
the  randomized  decision  rule.  The  decision  rule  g(2)  can  be 
summarized  as  follows 
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q{^)  ^  Prob.{Di|jr}  (2.8) 

I  1  if  1(2)  >  c 
_  ;  8  if  Jl(x)  ■  c 
t  0  if  S.{^)  <  c 

The  operating  values  of  c  and  8  are  chosen  to  satisfy  the 
requirements  of  the  detection  criterion  being  used.  If 
Prob.  {Jl(y)«c}  »  0  the  value  of  8  is  of  no  significance. 
From  now  on  we  will  not  consider  randomized  decision  rules 
and  will  assume  8  ”  1*  Any  monotone  increasing  function  of 
the  likelihood  ratio  is  an  equivalent  detection  statistic. 
We  will  frequently  use  the  log-likelihood  ratio  z(y)  defined 
below 


z(y)  ^  ln(a(y))  (2,9) 

2.3  Simple  and  Composite  Signal  Hypotheses 

When  both  ft(x)  and  fi>(y)  are  completely  known  the 
hypotheses,  Ho  and  Hi,  are  termed  simple.  In  this  case  the 
optimum  detection  statistic  £(y)  or  z(y)  may  readily  be 
computed.  However  when  the  signal  has  some  unknown 
parameters  e,  the  Hi  hypothesis  is  termed  composite 
[15,  p.  86].  If  the  noise  has  no  unknown  parameters  the  Ho 
hypothesis  remains  simple.  We  will  be  interseted  in  the 
situation  where  the  H|  hypothesis  is  composite  and  the  Ho 
hypothesis  is  simple.  When  H|  is  composite  it  is  assumed 
that  the  observer  knows  the  conditional  p.d.f.  A 

Bayesian  observer  based  on  the  state  of  his  knowledge  about 
e,  assigns  e  a  prior  p.d.f.  f(e).  The  marginal  p.d.f.  may 
now  be  computed 

■  Sq  fi(lll)£(i)<iO 


(2.10) 
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This  reduces  the  composite  signal  hypothesis  problem  to  a 
simple  signal  hypothesis  problem.  The  optimum  detection 
statistic  £(£)  is  now  obtained  as  follows 


iil) 


f  i(z) 

W 


•  /e  A(i|e)f(i)di 


(2.11) 


where 


a(^le)  ^  f i{i|e)/fo(i)  (2.12) 

i(y|e,)  is  the  conditional  likelihood  ratio.  iCjr)  in  this 
case  is  the  average  of  the  conditional  likelihood  ratios. 
Eq.  2.11  may  be  obtained  in  a  different  way,  fi(£il)  can  be 
written  as 


fi()tli)  ■  f i(2)f (ilx)/^(i) 


(2.13) 


Substituting  eq.  2.13  in  eq.  2.12  gives 


a  (III) 


f  i(z)f(llX) 


(2.Ua) 


or 

A(lll)^(i)  •  A(x)*(lll)  (2.14b) 

Rewriting  eq.  2.14b  gives 

A(l)  ■  A(xli)Hi)/f(ilx)  (2.15) 

Eq.  2.15  is  the  statement  of  the  Bayes^Birdsall  theorem. 
Integrating  both  sides  of  eq.  2.14b  with  respect  to  e  gives 


(2.16) 
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In  certain  situations,  even  for  a  Bayesian  observer  it  might 
be  unrealistic  to  assign  a  distribution  to  the  unknown 
parameter  vector  In  such  situations  an  approach  that  is 
often  used  is  to  fix  some  detector  and  then  analyze  the 
performance  of  the  detector  for  various  values  of  0.  The 
detector  may  or  may  not  be  based  on  the  knowledge  of  the 
actual  value  of  0.  In  Chapter  IV  we  will  use  this  approach 
extensively. 

2.4  Performance  Evaluation 


Given  the  decision  rule  g(y},  the  probabilities  of 
detection  P^  and  false  alarm  Pp^  can  be  determined  as 
follows 


Pjj  •  E(g{i)|Hi)  ■  Prob.{£(y)  >  c|Hi) 

•  ,  g(y)fi(2)di 

4(y)>c 


(2.17) 


/*  f(AlHj)d2 


and 


Pp^  *  E(g(y)jMo)  ■  Prob.{2(2)  >  clHo) 

-  ,  g(l)fo(l)ai 

ft(y)>c 

•  /,  ,  «o(iL)di 

B(y)>c 


(2.16) 


/*  f(A|Ho)d2 
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where  f(£|H^)  is  the  conditional  p.d.f.  of  the  likelihood 
ratio  under  hypothesis  K^.  Since  given  a  hypothesis  the 
decisions  D©  and  Dj  are  mutually  exclusive  and  collectively 
exhaustive,  we  have  the  following 

Pq  +  Pj,  »  1  (2.19a) 

It  follows  that  Pq  and  Pp^  are  sufficient  to  characterize 
receiver  performance.  From  eqs.  2.17  &  2.18  we  see  that  P^ 
and  Pp^  are  parameterized  by  the  single  parameter  c.  A  plot 
of  vs.  ^FA  for  all  possible  operating  values  of  c  is 
called  ihe  Receiver  Operating  Characteristic  (R.O.C.) 
curve.  The  R.O.C.  curves  completely  determine  performance 
for  all  operating  values  of  c.  R.O.C.  curves  have  been 
extensively  studied  by  Birdsall  [183. 

2.5  Two  Examples 

We  briefly  consider  two  examples  of  detection.  The 
first,  known  signal  in  Gaussian  noise,  introduces  the 
important  family  of  normal  R.O.C.  curves  and  also 
introduces  the  concept  of  normal  detectabilty .  The  second, 
Gaussian  signal  in  Gaussian  noise,  is  important  because  it 
occurs  frequently  in  the  later  chapters. 

2.5.1  Known  Scalar  Signal  in  Gaussian  Noise 

The  observation  under  Ho  and  Hj  is 

Ho  :  y  -  n  -  N(0,ci^  (2.20a) 

Hi 


:  y  »  n  +  s  ,  s  a  known  constant 


(2.20b) 
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Both  hypotheses  are  simple;  the  p.d.f,  of  y  under  Ho  and 
Hi  is 


f(ylHo)  =  (2iro^)~^'^^exp(-  y*/(2fl*)) 


(2.21a) 


f(y|Hi)  =  {2Trc^)”'^'^^exp{’-(y-s)*/(2a^)) 


(2.21b) 


the  likelihood  ratio  and  the  log-likelihood  ratio  are 


A(y)  =  exp((2ys  -  s^)/(2a^)) 


(2.22a) 


2(y)  »  (2ys  -  s*)/(2o*) 


(2.22b) 


We  define  d.  the  normal  detectability,  as  follows 


^  4  [E(z|Hi)  -  E(z|Ho)] 
vaFTTzJ 


(2.23) 


where  varo(z)  is  the  variance  of  z  under  Ho.  For  simple 
Gaussian  nypotheses  the  expression  for  normal  detectability 
simplifies  to 


d  ^  E(z(Hi)  -  E(z|Ho) 


(2.24) 


We  also  define  d*  as  follows 


d'  6  (2. 26) 

Physically  d  is  related  to  the  output  signal  to  noise  ratio 
(SNR)  and  is  used  as  a  performance  index.  For  a  fixed  value 

of  Pp^  larger  values  of  d  result  in  larger  values  of  Pq. 

2  ,  2 

For  the  above  problem  d  «  s  /o  and  the  statistics  of  z 
under  Ho  and  Hi  are 
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Ho  :  z  “  N(-d/2,d)  (2.26a) 

Hi  s  z  “  N(d/2,d)  (2.26b) 

The  probabilities  of  detection  and  false  alarm  are  given  by 
pQ  »  «(X+d’)  (2.27a) 

and 

Pp^  -  «(X)  (2.27b) 

where 

$(x)  -  (2ir)“^/^  exp(-t*/2)dt  (2.28) 

and  X  is  a  function  of  the  threshold.  The  R.O.C.  curves  so 
obtained  are  called  the  normal  R.O.C.  curves.  They  are 
usually  plotted  on  normal-normal  probability  paper  where 
they  plot  as  straight  lines  with  a  slope  of  1.  Normal  R.O.C 
curves  are  obtained  when  the  hypotheses  are  simple  and  the 
log-likelihood  ratio  has  a  normal  distribution  under  both 
hypotheses*  For  comparison  purposes  other  R.O.C.  curves 

are  also  plotted  on  normal-normal  probability  paper.  In 
fig.  2.2  a  family  of  normal  R.O.C.  curves  is  plotted  for 
several  d*  values.  In  later  chapters  we  will  use  d  and  d* 
to  compare  the  performance  of  various  detectors. 

2.5.2  Gaussian  Signal  in  Gaussian  Noise 

Now  assume  that  the  observation  is  a  vector  y  of  length 
L 

Ho  :  y  “  n  (2.29a) 

Hi  yv-n  +  s.nlis  independent  (2.29b) 


-r,  ■fv .iv  <•.  <•.  v.  < 
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where  hot  .  n  and  s  are  zero  mean  Gaussian  random  vectors  of 
length  L.  Let  and  Rg  denote  the  noise  and  signal 

autocorrelation  matrices  respectively 

R^  =  E(n.n'^)  (2.30a) 

and 

R_  »  eCs.s*^)  (2.30b) 

5  — »  • 

The  density  functions  of  the  observation  under  the  two 
hypotheses  are 

f(y|Ho)  »  (2ir)“^/2|R^r^/2exp(-/R;;^y/2)  (2.31a) 

and 

f(2|Hi)  ■  (2iT)“^/^(Rjj+Rgr^/^exp(-y'^(Rj^+RQ)"''y/2)  (2.31b) 

Where  |R|  >  determinant  of  R.  We  have  assumed  that  the 
appropriate  inverses  exist.  The  likelihood  ratio  and  the 

log-likelihood  ratio  are 

jl(y)  -  |Rnl^^^|Kn+R8r^'^^exp(^'^(R;;’-(Rj^+Rg)"bx/2)  (2.32a) 

z(y)  -  ^  z'^(R“’-(Rn+Rs)’bz  -  j  ln(  |R^+Rgl/|R„1 )  (2.32b) 

Since  the  second  term  in  eq.  2.32b  does  not  depend  on  the 
reception  an  equivalent  detection  statistic  is 

h(l)  -  Z^'tRn^  - 

the  above  may  be  written  as  follows 
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Ti(i)  =  (2.34) 

but  linear  minimum  mean  squared  error 

(m.m.s.e.)  estimate  [15]  o£  the  signal  vector;,  denote  this 

A 

by  s,  then 

T\iz)  ®  (2.35) 

This  result  will  be  used  extensively  in  Chapter  V. 

A  commonly  used  sub-optimum  detector  for  the  Gaussian 
signal  in  Gaussian  noise  case  is  the  Energy  Detector.  The 
energy  detector  uses  the  total  energy  in  the  observation  as 
the  detection  statistic.  Let  e(;^)  denote  the  detection 
statistic  for  the  energy  detector  then 

e(2)  »  (2.36) 

Note  that  the  optimum  detector  simplifies  to  the  energy 
detector  when  both  signal  and  noise  are  white.  In  Chapter  V 
we  will  often  compare  the  performance  of  the  optimum 
detector  with  that  of  the  energy  detector. 


CHAPTER  III 


PROBLEM  FORMULATION  I 

In  this  chapter  we  first  describe  the  underwater  | 

acoustic  channel  and  the  acoustic  propagation  model.  This  | 

is  followed  by  a  description  of  the  receiver  geometry.  | 

Finally  we  discuss  the  analytic  development  of  the  detection  | 

problem.  The  analytic  development  is  carried  far  enough  to  | 

allow  the  formulation  of  the  abstract  problems  discussed  in  | 

Chapters  IV  and  V-.  | 

3. 1  The  Underwater  Acoustic  Channel  j 

Models  similar  to  the  model  developed  here  for  the  | 

underwater  acoustic  channel  have  been  used  by  several  | 

authors  [5,10,19,21].  The  model  incorporates  a  uniformly  N 

moving  acoustic  source,  a  fixed  receiver  and  multipath  | 

propagation.  The  nth  path  in  the  channel  is  assumed  to  i 

exhibit  a  real  or  complex  gain  c^,  where  the  path  gain  ] 

magnitudes,  |c^|,  are  assumed  to  be  normalized.  We  assume  a  \ 

slowly  varying  ocean  so  that  the  c^^'s  and  the  number  of  n 

paths,  N,  with  significant  gain  remain  essentially  constant  t 

over  the  observation  interval.  The  propagation  delay  of  the  I 

nth  path  is  denoted  T^it) ,  where  j 

! 

’n<‘>  ■  ■'n  * 

I 


T-  is  the  path  delay  at  t«0,  either  the  start  or  the  middle 


‘n 

of 


the  observation  interval. 


is  the  path  delay 


derivative  or  the  doppler  shift  ratio.  For  receding  sources 


T^>0  and  for  approaching  sources  tJ^<0.  is  of  the  order 
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of  v/c  where  v  is  the  radial  component  of  the  source 

velocity  and  c  is  the  speed  of  sound.  The  spread  in  and 

t'  is  of  the  order  of  0.5%  to  1%  of  max  and  max  t' 

n  n 

respectively  [10j.  We  can  view  eq  3.1  as  the  truncated 
Taylor  series  expansion  of  Tj^(t).  This  models  uniform 
source  motion,  to  model  more  complicated  source  motion 
eq  3.1  would  have  to  be  augmented  by  including  higher  order 
terms  from  the  Taylor  series.  Putting  the  above  information 
together  the  underwater  acoustic  channel  can  be  viewed  as  a 
linear  filter  with  impulse  response  h(t),  where 


h{t) 


N 

n»1 


Tn(t)) 


(3.2) 


■  II 

n«i 


where 


'  -■‘k 


(3.3) 


According  to  the  above  model,  if  a  source  transmits  a  signal 
s(t),  the  receiver  will  receive  a  signal  r(t)  given  by 


N 


‘  ■  ''n> 


(3.3) 


To  complete  the  description  of  the  acoustic  channel  we 
incorporate  additive  noise  n(t)  in  the  reception,  so  that 


r{t) 


N 


(3.4) 


.>  ».V  «  I,-.  V  «  .V  y.  .  •, 
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We  assume  that  the  noise  is  independent  of  the  signal*  We 
also  assume  that  the  noise  is  a  zero  mean,  stationary, 
Gaussian  process.  The  Gaussian  assumption  is  widely  used  in 
the  underwater  acoustics  literature  [2,6,11]. 

3.2  Receiver  Geometry  and  Hypothesis  Formulation 

The  receiver  geometry  used  is  given  in  fig.  3.1.  Two 
receivers  RCVR1  and  RCVR2  are  placed  a  distance  d  apart.  In 
a  more  general  setting,  where  the  interest  is  in  localizing 
the  acoustic  source  in  three  dimensions  and  estimating  its 
velocity,  more  than  two  receivers  are  used.  We  will  assume 
that  the  distance  between  the  acoustic  source  and  the 
receivers  is  of  the  order  of  150  kilometers.  Assuming  that 
the  speed  of  sound  in  the  ocean  is  of  the  order  of  1500 
meters/second,  the  propagation  delay  from  the  source  to  the 
receiver  will  be  of  the  order  of  100  seconds  and  the 
differential  path  delays  will  be  of  the  order  of  tenths  of  a 
second.  We  define  T^  as  the  mean  travel  time  of  sound 
between  the  two  receivers 

T^  A  d/c  (3.5) 

where  c  is  the  mean  group  speed  of  sound  in  the  ocean.  The 
difference  in  mean  time  delays  to  RCVR1  and  RCVH2  is 
approximately  T^. 

We  assume  that  both  the  receivers  have  a  finite 
observation  interval  of  T  seconds.  We  will  assume  that 
baseband  signals  are  available  to  both  the  receivers.  The 
complex  demodulation  that  might  be  necessary  to  obtain  the 
baseband  signals  will  not  be  discussed.  The  receptions  at 
RCVR1  and  RCVR2  under  Ho  are 


Ho  :  yi{t)  -  ni(t)  ,  -T/2<t<T/2 


(3.6a) 
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Fig*  3.1  Receiver  geometry. 


Ho  :  yalt)  ■  n2(t)  ,  -T/2<t<T/2  (3.6b> 

the  n^(t)  are  stationary,  complex,  Gaussian  and  zero'-mean 
processes.  We  assume  that  the  receiver  separation  d  is 
large  enough  so  that  ni(t)  and  n2(t}  are  independent,  we 
also  assume  that  ni(t)  and  n2(t)  are  identically 
distributed.  The  noise  power  in  general  is  not  known. 
Under  the  signal-t^noise  hypothesis.  Hi,  we  have  the  following 
situation 


N 

Hi  :  yi(t)  =  ni(t)  +  ^Ci^s(aij^t  -  tij^)  ,-T/2<t<T/2  (3.7a) 

n=1 


M 

Hi  :  y2(t)  =  njCt)  +  ^C2jj,s(a2mt  -  12^^)  »-T/2<t<T/2 

m=i 


(3.7b) 


N  and  M  are  the  number  of  paths  to  RCVR1  and  RCVR2 
respectively.  Ci^^  and  C2jj  are  the  complex  path  gain 
coefficients,  ai^  and  ctz^  are  the  doppler  parameters  (as 
defined  in  eq.  3.3).  The  signal  s(t)  is  a  zero  mean, 
complex  and  stationary  Gaussian  random  process  on  the 
infinite  time  interval.  The  autocorrelation  of  s(t)  is 
denoted  Rg(T),  where 

Rg(T)  ■  aR(T)  (3.8a) 


and 


1R(t)|  <  1  (3.Bb) 

We  assume  R(t)  is  known,  however  "a"  in  general  is  unknown. 
We  also  assume  that  the  stochastic  signal  and  noise  are 
independent  and  are  individually  and  jointly  ergodic  [22], 
The  ergodic  assumption  allows  us  to  approximate  sample 
averages  by  long  time  averages  and  vice-versa.  We  assume 
the  signal  process  is  broadband,  so  that  the  signal 
autocorrelation  is  relatively  narrow  compared  to  the 
differential  path  delays.  For  our  purposes  a  signal  with  a 
bandwidth  of  the  order  of  10  Hz  or  more  is  considered 
broadband. 

3.3  Post  Reception  Processing  for  Detection 


The  usefulness  of  a  correlator  detector  has  already 
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been  pointed  out  in  the  introduction.  To  reiterate;  the 
correlation  allows  us  the  freedom  to  introduce  independent 
delay  and  doppler  search  parameters.  In  our  study  we  seek 
to  match  the  differential  dopplers  and  the  differential  path 
delays.  As  pointed  out  in  the  introduction  we  will  use  the 
Bivariate  Normalized  Crosscorrelation  (BNC)  functions 
Yip(3fT)  as  the  observable  for  the  detector. 


1  • 

/  yi($t)y2(t  -  T)dt 


“T/2 


(3.9) 


1  T/2  2  1 /o  1  •  21 /o 

(?  /  ,  |yi(t)|  ;  |y,(t)|  at)’/* 

A  -T/2  ‘  -T/2 


A  plot  of  the  magnitude  of  the  BNC  function  for  various 
values  of  6  and  x  is  termed  the  ambiguity  diagram.  Peaks 
are  obtained  in  the  ambiguity  diagram  when  the  search  values 
of  6  and  x  match  the  actual  values  of  differential  doppler 
and  delay  for  a  path  pair.  Our  goal  is  to  combine  the 
information  contained  in  the  peaks  in  a  manner  that  allows 
us  to  make  inferences  about  the  presence  or  absence  of  the 
source.  A  block  diagram  of  the  post  reception  processing  is 
given  in  fig.  3.2. 

We  now  look  at  the  BNC  function  in  more  detail*  As  in 
[10]  we  assume  that  the  integration  or  observation  times  T 
are  long  enough  for  the  denominator  of  to  stabilize. 
That  is,  we  assume  T  is  long  enough  for  the  individual  power 
measurements  to  stabilize.  Once  the  power  measurements 


stabilize  it  is  enough  to  consider  the  numerator  of  Yip(0*T)» 
denoted  by  Yy^  (B»t)»  as  being  representative  of  Yt(0»y). 

We  first  examine  the  denominator  terms  of  y<i>(0#y). 
Under  Hq  we  obtain 


i  I 

?  if/2 


lYilt) 


1 


|n^(t) 


(3.10) 
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Pig.  3.2  Block  diagrara  of  the  post  reception  processing 
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using  the  ergodic  assumption  we  obtain 


1  T/2  2 

^  ;  lni(t)|  dt  s  R  (0)  ,  T  large 
i  -T/2  ^  "x 

Under  Hi  we  obtain 


(3,11) 


1  t/2  2  1  2 

f  ,  lYiCt)!  dt  -  ^  /  |n.(t)l  dt 

i  -ip/O  *  *  -T/9  * 


T/2 


-T/2 


T/2 


,  t/2  g 

2  L-'in  f  /j/j  ni(t)s(oj„t  -  Ti„)dt 


l=<“in‘  -  ’in'l 


•t/2 


dt 


(3.12) 


*  J_  ’  lr/2 

Since  the  source  speeds  are  very  small  compared  to  the  speed 
of  sound,  it  is  reasonable  to  approximate  the  o^j  by  unity 
'in  the  denominator  integrals.  This  allows  us  to  use  the 
ergodic  assumption  in  eq.  3.12.  For  large  T  consider  the 
terms  of  eg.  3.12  individually 


1)  i  lni(t)|*dt  «  R„  (0) 
^  -T/2  ^  "i 


(3.13a) 


2)  2  Re  )  ^  /  nJ(t)8(aj„t  -  tj,^)dt  »  0  (3,13b) 

n  -T/2 
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V — .  2  1  T/2 


(3.13c) 


n 


V—  *  1  'r/2 

4)/  /  c>  c*  s!  r  s(a*  t  “ 

'/  /  in  im  T  “'“in'- 


n  m/n 


-T/2 


T*  )s*(a-  t  -  T-  )dt  (3.13d) 
^in'  ‘im 


-  ^  ^  c«  c*  R  (t*  ”  t*  ) 

_  /  in  im  s  im  ‘in' 


n  m^n 

From  the  broadband  assumption  we  have  min  jtjjjj  -  “'^inl^'^s' 

m/n 

;#here  is  the  signal  correlation  time,  so  that 

^s^'^im  **  ^in^“^*  given  the  different  phases  for  the 

c^j's  the  terms  in  eq.  3.13d  add  up  incoherently.  Given  the 
above  it  is  reasonable  to  use  the  approximation 


H  II=in'=tm‘'s<nin  '  Un>  “  “ 

n  m?n 


(3.14) 


Putting  the  approximations  of  egs.  3.13  and  3.14  together  we 
obtain 


1  T/2  2  r —  2 

Hi  !  f  /  |yi(t)|  at  =  R„,(o)  +  )  Icinl  *s<“> 

*T/  it  X 


For  the  approximation  to  be  valid  we  require  the  integration 
time  T  to  be  much  larger  than  max(T|^,Tg),  where  and  Tg 
are  the  noise  and  signal  correlation  times  respectively.  In 
general  Tg>Tj^,  so  we  require  T  »  Tg.  For  large  but  finite 
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T  the  right  hand  side  of  eg.  3.15  is  the  mean  value  of  the 
integral  on  the  left  hand  side.  We  expect  a  fluctuation 
about  the  mean,  we  assume  that  T  is  long  enough  so  that  the 
standard  deviation  is  much  smaller  than  the  mean. 

We  now  consider  the  numerator  of  the  BMC  function, 
Y,,n  (Bfi).  In  general  Y„n  U»t)  is  a  function  of  the 

Uuip 

integration  time  T,  the  time  delay  correction  t,  and  the 
doppler  correction  3. 


A  1 

Y,,n  ^  ,  yi(pt)yi(t  -  T)dt 


-T/2 


(3.16) 


Under  Ho  we  obtain 


1  * 

Ho  i  Y,,n  "  !?  /  ni(3t)n2(t  -  T)dt 


-T/2 


(3.17a) 


since  ni(t)  and  n2(t)  are  independent  and  zero  mean  we 
obtain 


Ho  :  Ynn  ^Brx)  «  0 


(3.17b) 


Under  H^  eg.  3.16  has  four  terms,  listed  below  individually 


1  * 

Dm/  ni(3t)n2(t  -  T)dt 

*  -T/2 


T2m  “  T)dt 


(3.18a) 


(3.18b) 


r —  .  T/2  * 

D>n  ?  L/,  ■ 


T)s(S«ij^t 


Ti„)dt 


(3.18c) 


. 


V.  W  V 
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(3.18d) 


.s*(a2^t  -  T2m  -  T)dt 

Long  time  integration,  T  »  max(T„,T_),  results  in 
egs.  3.18a,  3.18b  and  3.18c  approaching  zero  in  the  mean 
squared  sense.  All  three  terms  (3.18a,  3.18b  &  3.18c)  may 
be  considered  as  additive  "noise"  in  the  observation, 
although  the  "noise"  is  no  longer  independent  of  the 
"signal."  We  need  to  examine  the  behavior  of  eq.  3.18d 
under  long  time  integration.  Consider  the  following  value 
of  8: 
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^  T/2 


T,„)S  (oj„t 


4* 

When  T  =  Ti£  -  Taj  the  ci.CjjRgCx  -  (tij  -  peaks. 
Similarly  we  can  define  a  total  of  MN  $  values  for  which 
expression  like  eg.  3.20  are  obtained.  Since  we  know  the 
functional  form  of  t  terms  similar  to  the  first  term  in 
eg.  3.20  can  be  considered  the  ’’signal  term”  and  the  rest 
can  be  considered  a  part  of  the  additive  noise.  Here  we 
again  use  the  broadband  assximption,  that  is,  the  range  of  t 
for  which  Rg(T  -  (tij  ”  Taj))  is  significant,  the  other 
terms  in  eg.  3.20  are  small.  The  problem  then  becomes 
similar  to  the  set  of  problems  called  "signal  known  exactly 
except  for  specified  parameters  plus  additive  noise"  [16]. 
The  observation  is  two  dimensional  in  that  there  are  two 
search  parameters,  B  &  t.  This  is  a  non-trivial  problem  to 
solve,  to  begin  with  we  would  have  to  determine  the 
distribution  function  of  the  terms  lumped  together  as  noise. 

Work  has  been  done  by  several  authors  [23,24,25,26]  in 
the  area  of  determining  the  density  function  of  the  output 
of  an  analog  crosscorrelator  when  the  inputs  are  stochastic. 
However  our  formulation  is  complicated  in  that  the  received 
signal  is  the  sum  of  several  delayed  and  doppler  shifted 
versions  of  the  source  signal.  At  this  stage  the  custom  is 
to  remove  the  multipath  assumption  and  solve  a  problem  based 
on  a  single-path  assumption.  We  choose  to  keep  the 
multipath  assumption  intact  and  select  solvable  abstract 
problems,  derived  from  the  above  problem,  for  further  study. 
The  results  obtained  from  the  study  of  the  abstract  problems 
will  allow  us  to  form  rules  of  thumb  and  guidelines  for  an 
experimenter  faced  with  the  general  problem  formulated  in 
this  chapter. 


CHAPTER  IV 


EXTENDED  M’ORTHOGONAL  SIGNALS  AND  THE  FIRST  ABSTRACTION 

In  Chapter  III  we  saw  that  a  solution  to  the  exact 
problem  as  formulated  was  non-trivial  and  perhaps  not 
possible.  In  this  chapter  we  study  the  first  of  two 
abstractions  of  the  problem  formulated  in  Chapter  III.  The 
object  of  the  abstraction  in  this  chapter  is  to  find  out  how 
many  of  the  unknown  3,t  parameter  values  have  to  be  matched 
for  good  or  acceptable  performance. 

4. 1  Abstraction  Formulation 

In  practise  the  ambiguity  surface  is  usually  quantized; 
being  evaluated  at  discrete  values  of  3  and  t.  The 
quantization  is  usually  "fine”,  that  is  the  steps  in  t  and  3 
are  chosen  small  enough  to  yield  a  continuous  appearing 
surface,  but  large  enough  to  make  the  computing  job 
possible.  For  the  purposes  of  analysis  we  will  make  the 
usual  simplification  that  the  quantization  is  coarser  than 
used  in  practise  and  that  the  quantizations  are  "matched"  in 
the  sense  that  the  response  of  a  path-pair  falls  on  one  and 
only  one  point,  and  the  value  there  is  due  to  perfect 
alignment  in  time  and  frequency.  We  call  the  grid  point  a 
"cell"  and  the  computed  ambiguity  value  the  "cell  value". 
For  simplicity  the  cells  are  integer  indexed,  j«1,2,...,L 
and  the  cell  values  aie  denoted  by  yj.  The  total 
observation  is  thus  a  list  or  vector  l"y»,...»y£^  of  length 
L. 

The  observation  model  will  be  intentionally 

oversimplified  at  the  cell  statistics  level.  Specifically 
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under  Ho,  the  cell  values  are  i.i.d.  N(0,1).  Under  Hi,  C 
of  the  cells  are  i.i.d.  N(s,1)  and  the  rest  of  the  cells 
are  i.i.d.  N(0,1).  C  is  the  number  of  path  pairs  and  in 
the  notation  of  Chapter  III  C»M.N. 

The  complications  arise  from  the  composite  nature  of 
Hi.  Let  e  be  the  set  of  cells  whose  mean  is  s  under  Hi;  e 
is  the  "signal  parameter",  a  set  of  cell  indices.  The 
number  of  possible  0  is  nQ»L!/C! (L-C) ! .  We  refer  to  0  as  a 
"pattern"  or  "signal"  and  as  the  number  of  possible 
patterns  or  signals. 

Following  the  Bayesian  philosophy  we  put  a  prior 
distribution  on  the  signal  parameter.  If  the  prior  is  piQ^) 
where  i»1 ,2, . . . ,ng,  then  the  likelihood  ratio  is 

"0 

Jl(y)  -  y~^exp(2(y|0.  ))p(0.)  (4.1) 
where 


z(y|0j)  ■  s.^~  yj  -  .SCs*  (4.2) 

370i 

note  z(y|0^)  has  a  normal  distribution  and  exp(z(yj0^))  has 
a  log-normal  distribution  [29]. 

The  first  difficulty  one  would  encounter  if  one  were  to 
try  to  base  a  receiver  on  jt(y)  as  given  in  eg.  4.1  is  the 
size  of  the  parameter  space.  For  a  well  localised  signal 
the  ambiguity  plane  might  be  modelled  by  as  small  as  a  100 
by  100  square  i.e.  L«10,000.  For  a  4  path  by  4  path 

propagation  the  value  of  C  is  16;  there  are  n^-IO®*^*^^ 

possible  signals  and  as  such  2(y|0j^)  have  to  be 

evaluated  before  the  likelihood  ratio  can  be  calculated. 
The  primary  difficulty  lies  in  the  enormity  of  the  number  of 
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possible  signals  coupled  with  the  fact  that  we  lack  a  simple 
structure  to  reduce  this  to  a  sequence  of  smaller  sub¬ 
hypotheses.  Of  course  we  do  not  want  to  build  such  a 
receiver,  we  only  want  to  evaluate  its  performance.  We 
immediately  encounter  the  second  type  Cyf  difficulty;  namely 
that  the  determination  of  the  distribution  of  Jl(y)  is 
seemingly  hopeless  because  the  z(y|0p  are  not  independent 
over  i. 

Since  we  have  no  simple  way  of  determining  the 
performance  of  the  optimum  receiver  wo  must  consider  an 
alternative  approach  that  will  allow  us  to  make  statements 
about  the  problem  posed  at  the  beginning  of  this  chapter, 
i.e.  how  many  of  the  8,i  parameter  values  have  to  be 
matched  for  good  or  acceptable  performance.  Now,  out  of  the 
n»  possible  signal  patterns,  we  could  select  at  most  L/C 
that  do  not  have  a  cell  in  common.  For  the  simple  numerical 
example  of  the  100  by  tOO  square  and  a  4  path  by  4  path 
propagation,  such  a  partition  of  the  ambiguity  plane 
contains  only  625  patterns.  It  is  worthwhile  to  note  that 

A  *1  AM  « 

there  are  more  than  su'jh  partitions. 

In  the  following  we  consider  a  receiver  based  on  such  a 
partition,  i.e.  based  on  orthogonal  signals.  This  is  done 
primarily  for  mathematical  convenience;  however,  it  is 
realistic  in  that  the  receiver  designer  is  also  faced  with 
the  enormity  of  n^  and  is  forced  to  simplify.  For  example  a 
receiver  designer  may  falsely  believe  that  the  doppler  is 
the  same  on  all  the  paths  to  a  given  receiver,  in  other 
words  the  designer  o^'erlooks  differential  doppler.  This 
means  that  the  source  will  affect  only  one  narrow  p  band  of 
the  ambiguity  plane  (in  our  simplified  ambiguity  grid  this 
band  corresponds  to  one  p  index,  that  is,  the  doppler  is 
centered  and  the  frequency  resolution  sharp  enough) .  To 
match  our  evaluation  model  exactly  the  designer  would  have 
to  know  all  of  the  differential  delays  well  enough  to 
specify  the  C  patterns  without  overlap.  This  is 
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unrealistic r  so  the  following  evaluation  is  optimistic  and 
can  be  considered  an  upper  bound  on  performance. 

This  chapter  will  consider  a  receiver  based  on 
orthogonal  signals  and  examine  its  performance  under  several 
conditions;  under  the  assumption  that  the  partition  was 
correctly  chosen  and  the  true  signal  is  one  of  the  patterns 
in  the  partition,  and  under  two  conditions  where  the  true 
signal  is  not  one  of  the  patterns  in  the  partition.  These 
latter  conditions  are  the  ones  of  interest,  for  they  may 
shed  some  light  on  how  well  one  must  match  the  true  pattern 
to  perform  well. 


4.2  The  Orthoqonal-Siqnal-Based  Receiver 

In  this  section  we  examine  the  orthogonal-signal-based 
receiver,  introduced  in  the  previous  section,  in  more 
detail.  Given  a  partition,  we  first  reconfigure  the 
observation  y  into  a  matrix  so  that  each  in  the  partition 
is  one  row  of  C  columns.  Without  loss  of  generality  we  can 
re-index  the  patterns  so  that  the  receiver  is  based  on  the 
first  L/C  patterns  (for  simplicity  in  analysis  we  will 
assume  that  R>L/C  is  an  integer).  C  is  a  mnemonic  for  both 
the  number  of  cells  in  the  signal  pattern  and  for  the  number 
of  columns  in  the  reconfigured  ambiguity  diagram.  R  is  the 
mnemonic  for  the  number  of  rows  in  that  reconfiguration.  We 
also  restrict  ourselves  to  the  simple  case,  where  the 
receiver  is  based  on  the  traditional  worst  case  assumption 
that  the  R  signals  are  equally  likely.  Formally  we  base  the 
receiver  on  the  following  prior  distribution 


pCe^) 


1/R  ,  iff  i-1,2,...,R 

0  ,  iff  i«R+t,...,nQ 


(4.3) 


with  the  reminder  that 


e^nOj  is  empty  for  1<i<j<R 


(4.4) 


The  true  signal  pattern  will  be  denoted  by  0j  or  simply  by 
J.  Also  for  the  rest  of  this  chapter  Hi  will  mean  Hi  «  . 


4.2.1  The  M-Orthoqonal  Signals  Evaluation 


Here  we  evaluate  the  above  receiver  under  the 
assumption  that  the  true  signal  is  one  of  the  patterns  in 
the  partition;  that  is,  J<R.  This  case  is  known  as  "signal 
one  of  M-Orthogonal  Signals"  or  simply  M-Orthogonal  Signals, 
abbreviated  as  MOS  [13,18,27,28].  M  is  numerically  equal  to 
R  here.  The  usual  but  optimistic  normal  detectability 
approximation  is 


«MOS  ■ 


(4.5) 


4.2.2  The  Extended  M-Orthoqpnal  Signals  Evaluation 


Here  we  evaluate  the  orthogonal-signal-based  receiver 
under  two  conditions  where  the  true  signal  is  not  one  of  the 
patterns  in  the  partition.  We  call  this  Extended  M- 
Orthogonal  Signals  (BMOS)  as  the  receiver  is  designed  for 
MOS,  but  the  true  signal  is  not  one  of  the  design  signals. 
The  true  signal  does  have  cells  in  common  with  the  design 
signals.  Let  the  number  of  cells  that  the  true  signal  has 
in  common  with  the  ith  design  signal  be  denoted  by  K^,  then 

■  card{0jn0^}  (4.6) 


where  card{  }  is  notation  for  the  cardinality  of  the  set. 

For  further  analysis  we  make  the  following  assumptions 
and  restrictions 


(4.7a) 
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C  =  2 


a 


(4.7b) 


b  >  2a 


(4.7c) 


R 


(4.7d) 


Later,  we  will  be  interested  in  the  situation  where  L*8192, 
C«64  and  R»i28  (this  corresponds  to  an  8  path  by  8  path 
propagation  model). 

.  Consider  positive  integers  R  and  C,  R>C,  and  any  set  k 


of  R  non-negative  integers 


R 

such  that  /  K.  » 

i»1 


denote  the  number  of  decompositions  of  C 
summands  without  regard  to  order  [35  pg.  825]. 
obtain 


C.  Let 

into  integer 
For  C“64  we 


#{<}  -  1,741,630  (4.8) 

The  receiver  design  was  based  on  R  orthogonal  signals, 
each  signal  affecting  C  cells.  For  each  signal  strength,  s, 
there  are  possible  conditional  R.O.C.  curves,  one  for 
each  set  k.  It  is  unrealistic  to  consider  a  distribution  on 
{<}  and  seek  an  average  R.O.C,  An  upper  bound  on  these 
R.O.C. '8  is  given  by  the  MOS  evaluation,  where  ic>(Ki , . . . ,Kj^) 
is  specified  by 


K.  -  {  ^  (4.9) 

^  to,  i-2,3,...,R 

A  lower  bound  on  these  R.O.C. 's  is  given  by  the  sparse 
condition 

K  «  I  1  »  i-1,2,...,C  (4.10) 

^  to,  i«C+1,..,,R 


To  shed  light  on  the  central  question  "how  well  does 
the  receiver  have  to  come  to  matching  the  actual  signal?" 
two  types  of  intersecting  sets  were  selected  and  performance 
evaluated  for  C=64  and  R=128.  The  first  type  which  we  call 
standard  EMOS  is  given  by  the  condition 


2^  ,  i=1,2,...,2^"^ 
0  ,  o.w. 


(4.11) 


The  second  type  which  we  call  modified  EMOS  is  given  by  the 
condition 


,  i=1 

,  i»2,3, . . . ,65-Ki  (4.12) 

,  o.w. 

performance  for  standard  EMOS  should  be  better 
than  the  performance  for  modified  EMOS.  We  are  also 
interested  in  determining  how  much  the  performances  differ. 
That  is,  how  much  does  it  help  performance  if  the  largest 
occurs  several  times? 


Intuitively 


4.3  Receiver  Performance  Evaluation 

In  the  following  we  will  discuss  performance  evaluation 
techniques  for  the  standard  EMOS  in  detail.  Later  when  we 
discuss  Monte-Carlo  methods  for  performance  evaluation  we 
will  also  consider  modified  EMOS.  Unless  otherwise 
specified  by  EMOS  we  will  mean  standard  EMOS. 

Let  S  denote  the  number  of  rows  that  intersect  the  true 
signal  6j  and  let  K  denote  the  value  of  in  the  S  rows. 
We  can  now  write  down  the  likelihood  ratio  for  the  EMOS 
evaluation 
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S  R 

Z(y)  «  exp(Zj^(yl0. ))  +  ^  exp(zQ(y  |  0. ) ) )  (4.13) 

i7T  ra+1 

where  the  subscripts  K  and  0  are  used  to  denote  the  number 
of  cells  the  ith  row  has  in  common  with  the  true  signal. 
Now 


T —  2 

Zj^(yl0i)  «  s.) _  yj  "  ‘SCs  ,  i»1,...,S 

and 


^ —  * 

!Q(y(0^)  ■  ”  . 5Cs  f 


(4.14a) 


(4.14b) 


We  are  interested  in  the  statistics  of 
Hi 


yj  under 


Ho  and 


Ho 


^  '  yj  ■  **  N(0»C)  i.i.d.  t  i*1»»»»fR 


(4.15a) 


sR  ♦  2^  ,  i»1 , . . .  ,S 
El  ,  i«S+1,...,R 


The  distribution  of  2Q(yl0i)  and  2jj(yl0j)  easily  follow  from 
the  above 


2o(y|0i)  “  N(“  §s’,Cs’) 


(4.16a) 
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Ho  :  Zj^(y|Gj)  “  N(-  |s*,Cs*)  (4.16b) 

Hi  :  ZK(y|8i)  ~  N({K  -  |)s*,Cs*)  (4.16c) 

Since  ZQ(y|0.)  and  Zj^(y|e.)  have  a  normal  distribution 
exp(zQ(y|e£))  and  exp(Z|^(y|  6^) )  have  a  log-normal 
distribution  and  £(y)  has  the  distribution  of  the  sum  of 
independent  log-normals.  Before  proceeding  further  we  state 
a  few  pertinent  facts  about  the  log-normal  distribution.  A 
detailed  discussion  and  further  facts  may  be  found  in  [29]. 
Let  random  variables  Z  and  Y  be  related  as  follows 

X  «  InY  (4.17) 

with 


X  -  N(v,o*)  (4.18) 

then  Y  is  said  to  be  log-normally  distributed.  Me  denote 
this  as  Y  **  ).The  density  function  of  Y,  f(y)  is 

;  (5^)^/2.(oy)'^exp(-  ^.(ilUC^)*)  ,  y>0 

f(y)  -  (4.19) 

to  .  y<0 

The  kth  moment  of  Y,  E(Y^)  is  given  by 

2  2 

E(Y*‘)  «  exp(kv  +  (4.20) 

and  the  variance  of  Y,  Oy  is 

2  2  2 

Oy  ■  exp(2v  0  ).(exp(o  )  -  1)  (4.21) 


O  O.*--  •w’  •-*  ^ *  •  *  *  •  •  O  •  •  •  '•  »  *  •  VV  iV  •  *  A-  \  *•  *  »•%  ‘T ,  •  .  •*  >  •*'.  <**.  ^1*  * 
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To  evaluate  the  receiver  performance  we  need  to  find 
the  distribution  of  the  sum  of  independent  log-normal  random 
variables.  Since  no  closed  form  expression  exists  for  the 
density  function  of  the  sum  of  independent  log-normal  random 
variables  we  must  resort  to  other  techniques,  i.e. 
numerical  or  series  approximation  methods.  We  will  discuss 
several  ways  for  evaluating  receiver  performance. 

4.3.1  Approximation  to  the  Detectability  Index 

The  normal  detectability  d  or  its  square  root  d*  are 
accurate  measures  of  receiver  performance  when  the  R.O.C. 
is  normal.  In  this  case  the  R.O.C.  is  clearly  not  normal, 
however  we  will  use  normal  detectability  as  an  approximate 
measure  of  performance  [16]. 


5 

E(a(y)lHi}  •  g(^E(£j^(yl6^)|Hi) 


ra-^i 


(4.22) 


Where  Aj^( .  )«exp(Zj^( . ) )  and  JIq(  .  )>»exp(zQ( . ) ) .  When  the 
R.O.C.  is  normal  we  have  [30] 


C(£|H|)  •  exp(d) 


(4.23a) 


E(a}Ho)  -  1 


(4.23b) 


So  assuming  £(y)  results  in  a  normal  R.O.C.  we  obtain  for 
eq  4.22 


*  |.exp(d^(K))  +  ^ 


(4.24) 


Where 


^EMOS  ^  'detect,  for 

(K)  ^  detect .  for 
eq.  4.24  simplifies 


Extended  M-Orthog.  signals 

one  row  with  K  signal  cells 
to 


exp(dgj^Qg)  -  1  =  |.(exp(dx(K))  -  1) 
We  will  consider  two  cases: 


1)  Both  di(K)  and  are  large. 

2)  Both  di(K)  and  ci^j^Qg  are  small. 

With  both  di(K)  and  dg^^^g  large  we  obtain 

exp(dgj^Qg)  s  g.exp{di^K)) 

or  taking  the  logarithm,  we  obtain 


(4.25a) 

(4.25b) 


(4.26) 


(4.27a) 


dgMos  ~  di(K)  -  ln(R/S)  (4.27b) 

When  both  and  di(K)  are  small  we  can  expand  exp(di(K)) 
and  Taylor  series  and  ignore  the  higher  order 
terms.  Thus  we  obtain 


‘^EMOS  “  (4.28) 

The  results  in  eg.  4.27b  and  eg.  4.28  are  similar  to  the 
results  obtained  for  the  M-Orthogonal  signals  case  [30]. 
Now  di(K)  may  easily  be  calculated  from  the  expression  for  d 
and  the  statistics  for  2j^(y|8^)  given  in  eq.  4,16 
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,v 


i 


V 


s 


'2 


cli(K)  ^  E(zj^(y|ei)  |Hi)  ■  B(zK^y|®iHHo) 


(K  -  §)s' 


{-  §)S^ 


(4.29) 


=  Ks 


So  for  large  d  we  have 


•^EMOS  “  "  ^"1 


(4.30) 


We  are  interested  in  the  situation  where  L«8192,  C»64 
and  R-128.  For  the  MOS  evaluation  (KsC)  we  desire  that  the 
detection  be  almost  sure;  we  will  use  this  to  dictate  our 
choice  for  s*  From  experience  a  d'  e  8  results  in  high 
detection  probabilities  for  very  small  false  alarm 
probabilities  (also  see  fig.  2.2).  Substituting  in  eq.  4.30 
we  obtain 


64  7.1n2 

- ^ - 


or 


s  s  1.03 


(4.31) 


In  the  following,  for  convenience  we  let 
above  choice  of  variables  we  can  find  d 
These  values  are  tabulated  in  table  4.1. 


EMOS 


s«1.  With  the 
for  several  K. 


4.3.2  Formal  Evaluation  of  the  p.d.f .  of  (v) 


Our  interest  is  in  deriving  the  density  function  for 


£(y)  (eg.  4.13)  under  both  the  signal-t-noise  hypothesis  and 
the  noise  alone  hypothesis.  Under  the  noise  alone 
hypothesis  eq.  4.13  simplifies  to 


it 


1n->.V>.V. 


K 

di(K) 

‘^EMOS 

^EMOS 

D 

1 

.5 

.707 

2 

in 

.707 

1.92 

1.39 

8 

8 

5.23 

2.29 

16 

16 

12.53 

3.54 

32 

32 

27.84 

5.28 

64 

64 

59.15 

7.69 

Table  4.1  Table  of  normal  detectability  for  several 
K  values,  for  signal  power  s»1. 


R 

K  ^exp(zQ(y|9^)) 


■  g.exp(Zjj{yl0^) ) 


;^.f-(lnRa)  ,  a  >  0 

U  Su 


0  ,  o.v. 


(2irCa  )“^/^.exp(-  ■5^(lna“(K  -  |  -  InR))*)  ,  a>0 


f...  (a)  = 


0  ,  o.w. 


(4.35) 


for  C=64  and  R*128  the  above  reduces  to 


(5l)'/2.gl.e.p(-  -  (K-36.852))^,  _ 


(a)  * 


0  ,  ow. 


(4.36) 


Since  the  are  independent  and  identically  distributed 

for  fixed  K,  the  density  function  for  the  likelihood  ratio 
Jl(y)  is  the  convolution  of  the  density  functions  of  the 
w. 's.  Or  if  we  let  4>„  (v)  be  the  characteristic  function  of 

f„  (a)  and  (b.  (v)  be  the  characteristic  function  of  f.  (a) 
''k  ^k  ■‘•k 

we  have 


^w  ^w  jva)da 

k  k 


♦l  (V)  •  (♦,  )*.(♦«  (v))”'* 


fi  (a) 
^k 


H  -^-oo  tv).exp(-jva)dv 
0  ,  o.w. 


(4.37) 


(4.38) 


(4.39) 


(4.40) 


Iff  J*-®  (v).exp(-jva)dv 


fi  (a)  •* 


0  ,  O.W. 


(4.41) 


Now  if  we  choose  X  as  the  threshold  value  for  the  purposes 


cnMp 
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of  decision  making,  the  probabilities  of  detection  Pq(x|K} 
and  false  alarm  may  be  calculated  as  follows 

Pq(X|K)  =  /"  f;^^(a)da 

®  1  -  /q  fl  (a)da  (4,42) 

)c 

PpA  '  fx 

■  1  “  /q  f^^(a)da  (4.43) 


Formally  the  problem  is  solved.  However  to  get  the  R,0,C, 
curves  we  need  to  do  more  because  no  closed  form  solutions 
to  the  above  integrals{4,37  to  4,43)  exist.  We  need  to  find 
suitable  numerical  or  Monte-Carlo  procedures  for  evaluating 
the  R,0,C,’s,  In  the  following,  three  procedures  that  were 
considered  for  evaluating  the  R.O.C,  curves  are  described. 


4,3,3  Discrete  Fourier  Transform  Approach 


The  DFT  approach  for  evaluating  the  density  function 

for  the  svun  of  log-normal  random  variables  is  fairly 

natural.  The  procedure  for  evaluating  the  density  function 

and  the  R.O.C.  is  briefly  discussed. 

First  to  calculate  the  characteristic  functions  (v) 

k 

we  need  to  find  a  value  of  a#  say  n  such  that  for  all  K 


(4.44) 


Where  e  is  chosen  in  accordance  with  how  much  error  can  be 
tolerated  in  the  calculation.  The  next  step  is  to  sample 


f„  (a)  and 


(a) 

"O 


in  the  region  (O,^]).  The  sampled 


•v:?' 


sequences  are  then  padded  with  a  sufficient  number  of  zeros 
to  avoid  the  wrap  around  effect  in  later  calculation.  This 
is  followed  by  evaluating  the  DFT  of  the  two  sequences.  The 
DFT  of  the  sampled  versions  of  f„  (a)  and  f„  (a)  are  raised 

to  the  Sth  power  and  the  (R-S)th  power  respectively.  The 
sequences  so  obtained  are  then  multiplied  and  the  inverse 
DFT  operation  is  performed;  this  results  in  a  sampled 
version  of  f-,  (a).  A  similar  procedure  yields  the  sampled 

version  of  f-.  (a)*  A  suitable  numerical  integration 
J-o 


A  suitable  numerical  integration 


procedure  (say  Simpson's  rule)  then  yields  the  probabilities 
of  detection  and  false  alarm. 

In  our  problem,  for  the  values  selected  for  L,  C,  R  and 

s,  this  procedure  breaks  down.  From  the  values  given  in 

2 

Table  4.2  for  E(W|^),  El(W|^),  maxf^  (a)  and  the  value  of  a  at 

k 

which  the  maximum  occurs  we  see  that  the  f„  (a)'s  are  not 

''k 

well  behaved.  A  suitable  common  sampling  period  can  not  be 
found.  It  is  also  not  appropriate  to  model  the  densities  as 
6  functions  due  to  their  long  tails  as  evidenced  by  the  huge 
variance. 


This  section  is  based  on,  and  is  an  extension  of  the 
techniques  developed  by  Nolte  and  Jaarsma  [28]  for 
evaluating  the  R.O.C.  curves  for  the  detection  of  one  of  M> 
Orthogonal  signals. 

We  will  only  describe  the  calculation  for  Pq(x|K).  The 
techniques  for  calculating  are  the  same  except  that  S 
and  K  are  both  set  to  zero. 

We  can  write  the  4.  (v)  the  characteristic  function  for 

the  density  function  of  the  likelihood  ratio  in  terms  of  its 
real  and  imaginary  parts 


1-1 

o 


M 

» 


4) 


Xi 

(d 

•H 

u 

(d 

> 


§ 

•o 

c  • 

(d 

b  VO 


(U  (S 
JS  eo 

4J  » 
VO 


O  •• 
00 
m  » 
o  ^ 

•H  ^ 
C4 

(0  •» 


4J  ■» 

Id  o 

4J  H 
CO  t>d 


CM 


4) 
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(4.45) 


Substituting  the  above  in  eg.  4.39  gives 


fi  (a) 
•^k 


^  /*a>  [4>R^v)+j4»j(v)][cos(va)+jsin(va)]dv  ,a>0 
k  k 


(4.46) 


o.w. 


Since  f-,  (a)  is  real  eg.  4.46  reduces  to 
■^k 


f  1  (a) 


/-OD  C4R(v)cos(va)+4j(v)sin(va)]dv  ,oi>0 
k  k 


(4.47) 


,  o.w. 


Also  since  £■,  (a)  is  the  convolution  of  densities  that  are 
•^k 


zero  for  a  <  0,  f.  (a)  must  also  be  zero  for  a  <  0  i.e.  for 
“  -^k  • 


a  >  0,  f,  (-a)  ■  0.  This  implies  that 


/-«  4>R(v)cos(va)dv  •  (v)sin(va)dv  ,  a>0 

k  k 


(4.48) 


So 


£■1  (a) 
^k 


^  ♦r  (v)cos(va)dv  ,  a  >  0 
k 


0  f  o.w. 


(4.49) 


and  since  (v)  is  even 


(a)  • 
^k 


f  /q  ♦|i^(v)cos(va)dv  ,  o  >  0 


(4.50) 


,  o.w, 


Now  substituting  eg.  4.50  in  eg.  4.42  we  obtain 
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.  Pq(X|K)  =  ^  "  i  /o  ^0  <v)cos(va)dvda  (4.51) 

)c 

Changing  the  order  of  integration  in  eg»  4.51  gives 

Pq(X|K)  »  ^  "  f  /o  /o  cos(va)dadv  (4.52a) 

Ic 

Iv 

«  1  -  ^  /o  <t»R  (v)?i-^-^^idv  (4.52c) 

ic 


For  the  time  being  assume  that  d>o  (v)  is  available,  then  the 

integral  in  eq.  4.52c  can  be  calculated  using  Simpson's  rule 

between  the  zeros  of  the  function  The  number  of 

points  chosen  between  the  zeros  depends  on  how  fast  (v) 

varies.  The  integration  is  carried  out  until  successive 

—7 

approximations  to  the  integral  differ  by  less  than  10  . 

Now  we  develop  the  procedure  for  evaluating  (v)  at 

k 

the  selected  v  values.  We  can  write  the  characteristic  • 
functions  (v),  ^  (v),  and  (v)  as  follows 

W  A  W*.  *1. 


♦,,^(v)  »  l4>y^(v)  l.exp(jeo(s)) 

(4.53a) 

♦w  *  ^♦w  <v)  1  .exp( je)j(s)) 

k  k 

(4.53b) 

^v)  ■  |4t,  (v)  1  ,exp(  je,  (s)) 

(4.53c) 

Where  U»,jj(v)|#  0o(v),  (v)j,  0|j(v), 

K 

and  |4).  (v)| 
^  k 

,  0l(v) 

are  the  magnitude  and  phase  of  (v), 

"0 

(v)  and 
''k 

respectively.  Substituting  the  above  in  eq.  4.36  gives 


(v)l.exp(je3^(v)  »  (v)l^.|4»^^(v)|^“®’ 

Ic  Ic 


(4.54) 


,exp(  j(S9jj(v)  +  (R-S)0o(v))) 


It  follows  that 


♦r  (v)  »  |4>j^  (v)  1cos(62^(v) 
K  Ic 


(v)  I*.  U^^(v)  p'’®.cos(S0j^(v)  +  (R“S)eo(v) )  (4.55) 

Ic 


We  need  to  calculate  the  magnitude  and  phase  of  ^  (v)  and 

k 


(v)  at  the  required  v  values.  We  will  describe  the 

W  A 


procedure  for  obtaining  ^  (v);  4  (v)  is  obtained  using  an 

k  ® 


identical  procedure  with  K»0.  Now  a  (v)  is  given  by 


(V)  .  (2.0''''^  i.exp(-  ^(Ino  -  (K-§-lnR))’) 


.exp( jva)da 


(4.56) 


Define 


H 


K  -  S  -  InR 


(4.57a) 


and 


A  C 


(4.57b) 


Substituting  the  above  in  eq.  4.56  gives 


.exp( jva)da 


(4.58) 


Using  the  substitution  6  >  Ina  gives 
4>  (v)  =  (27ra*)'''^^  exp(“ 

Iv 

.exp(jve®)d3  (4.59) 

Since  we  need  to  start  the  numerical  integration  at  a  finite 
lower  limit  a  suitable  lower  limit  which  results  in 
negligible  error  must  be  found.  Let  x  desired  lower 

limit  then 

^^^(v)  ■  (2ircj*)’^'^^  exp(“^(^-^-!^)*)  .exp(  jve®)dB 

+  (2ir<3*)’^'^^  /“  exp(-^(^-j-!^)*) ,exp( jve®)de  (4.60) 
Then  the  error  E  is 

E  -  (2iro*)“^/^  exp(-  ^{i-5-!ii^)*).exp(jve^)dB  (4.61) 


and 


|E|  <  exp(-  (4.62) 

if  we  let  X  ■  ”6o  +  Uj{/  obtain 

|E|  <  (2iio')"'/^  exp(- 


(4.63a) 
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using  the  substitution  x  «  ^  we  get 

a 


|El  <  /I®  exp(-  x*/2)dx  e  10“* 


(4.63b) 


This  error  is  smaller  than  the  quantization  error  introduced 

on  the  machine  that  was  used  for  the  calculation. 

Let  (v)  and  (v)  be  the  real  and  imaginary  parts 
''kR  '^kl 

of  (j>„(v),  from  eq.  4.59  with  x  substituted  for  we  obtain 
k 


exp(- 

.cos(ve^)d6  (4.64a) 

exp{- 

.sin(ve^)d6  (4.64b) 

Eqs.  4.64a  t  4.64b  were  evaluated  using  a  5  point  Gaussian 
integration  procedure  [35].  From  -6a**’U|(  to  the  first  zero 

of  sfn^ve^}  the  5  point  gauss-quadrature  procedure  was 
applied  over  steps  of  a/4.  Thereafter  the  integration  was 

performed  between  successive  zeros  of  a®®{ve®).  The 
integration  was  continued  until  the  difference  betweeen 
successive  approximations  to  the  integral  became  less  than 

10  .  The  R.O.C.  curves  so  obtained  are  plotted  in 

fig*  4.1* 


4*3*5  Monte-Carlo  Methods  for  Evaluating  the  R.O.C* 


Monte  Carlo  simulation  methods  wore  used  to  confirm  the 


Pig.  4.1  R.O.C.  curves  Cor  the  Extended  M-Orthogonal 
signals  for  R*1,2,4,8r16>32. 
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fig.  4.2  R.O.C.  curves  for  the  modified  Extended 
K-Orthogonal  signals  for  R«1,2,4,8,16/12 
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results  obtained  in  fig,  4.1.  The  procedure  is  breifly 
discussed  here.  Basically  a  decision  making  device  is 
simulated.  The  device  sees  random  variables  with  the 
distribution  of  either  zCyjHo)  or  z(y|Hi).  According  to  a 
predefined  threshold  the  device  either  accepts  or  rejects 
the  Hi  hypothesis.  The  probability  of  false  alarm  for  a 
given  threshold  level  is  then  estimated  by  dividing  the 
number  of  z(y|Ho)  r.v.’s  for  which  Hi  is  accepted  by  the 
total  number  of  z(y|Ho)  r.v.’s  generated.  A  similar 
procedure  using  z(y|Hi)  r.v.'s  gives  the  probability  of 
detection.  The  procedure  is  repeated  for  as  many  threshold 
levels  as  points  are  desired  on  the  R.O.C.  curve.  Our 
simulation  evaluated  31  points  on  the  R.O.C.  curve  based  on 
10,000  trials  each. 

The  analysis  in  the  previous  was  restricted  to  what  we 
had  termed  standard  EMOS.  This  was  done  primarily  for 
analytic  and  computational  simplicity.  Monte-Carlo  methods 
allow  us  to  evaluate  the  R.O.C.  curves  for  the  modified 
EMOS,  This  was  done  for  the  following  parameter  values 
L=8192,  C=64,  R=128,  s=1  and  for  Ki=1 ,2,4,8, 16,32.  The 
R.O.C,  curves  so  obtained  are  given  in  fig,  4.2. 

4.4  Summary 

The  R.O.C,  curves  obtained  in  figs.  4.1  &  4.2  suggest 
that  when  the  design  signals  are  such  that  the  number  of 
signal  cells/row  that  are  common  to  the  true  signal  is 
small,  the  quality  of  detection  is  poor.  The  R.O.C.  curves 
in  fig.  4.2  in  some  sense  provide  a  lower  bound  to 
perfomance  in  that  we  assume  that  localizing  some  of  the 
signal  bearing  cells  does  not  provide  any  information  about 
the  other  signal  bearing  cells  in  the  ambiguity  plane.  We 
see  that  in  this  case  it  is  necessary  to  localize 
approximately  half  the  signal  bearing  cells  to  obtain 
performance  that  is  substantially  better  than  chance.  The 
R.O.C.  curves  in  fig.  4.1  assume  that  localized  signal 
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cells  provide  some  information  about  the  other  signal  cells 
in  the  ambiguity  plane.  We  see  that  performance 
substantially  better  than  chance  is  obtained  when  a  quarter 
to  a  half  of  the  signal  cells  have  been  localized.  The 
above  conclusions  are  based  on  the  fact  that  no  single 
signal  bearing  cell  dominates  the  ambiguity  diagram.  The 
EMOS  method  also  provides  us  with  a  technique  for  combining 
the  information  in  the  ambiguity  diagram  for  detection 
purposes. 


»  ^  ■»  _,  *  1.  •  —  *'1,  '  ^  *  •  *  •  •  '  •  •.  *«  *«  •_  *_  *  •_  ^  •*•*•*•*•**"»  »  « 


CHAPTER  V 


PERIODIC  RANDOM  SEQUENCES  AND  THE  SECOND  ABSTRACTION 

As  Stated  earlier,  the  MN  peaks  in  the  ambiguity 
diagram  are  generated  by  M+N  independent  sets  of  doppler  and 
delay  parameters.  Two  of  the  questions  addressed  in  this 
chapter  are 

1)  Should  the  receiver  treat  the  MN  peaks  as 
independent  entities  or  as  being  generated  from  M+N 
independent  parameter  sets?  In  other  words  should  the 
receiver  have  MN  degrees  of  freedom  or  M+N  degrees  of 
f  reedom? 

2}  Does  one  approach  offer  a  significant  improvement  in 
performance  over  the  other? 

There  is  no  simple  or  direct  way  of  answering  these 
questions.  Our  technique  is  to  study  a  solvable  abstraction 
that  captures  the  essence  of  the  above  problem.  To 
construct  the  abstraction  we  use  the  notion  of  Periodic 
Random  Sequences.  In  Sect.  5.1  Periodic  Random  Sequences 
[PRS]  are  introduced  and  a  framework  for  the  abstract 
problem  is  established.  The  following  sections  deal  with 
the  theory  of  detection  of  Periodic  Random  Sequences  in  some 
detail.  After  this  digression,  the  abstract  problem  is 
considerd  in  Sect.  5.4. 

5. 1  Periodic  Random  Sequences 

We  define  a  sequence  p(n)  to  be  periodic  random  with 
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period  P  if  it  satisfies  the  following  two  properties 


1.  p(n)  ~  N(0,o^)  i.i.d.  0  <  n  <  P-1  (5.1) 

2.  p(n)  >  p((n)modP)  all  n.  (5*2) 

Of  course  in  general  a  PRS  need  not  have  independent  or 
gaussian  samples  and  the  samples  may  not  necessarily  be  zero 
mean.  However,  unless  otherwise  stated,  by  a  PRS  we  will 
mean  a  sequence  that  satisfies  the  above  properties.  The 
autocorrelation  function  Rp(k)»E{p(n)p(n+k) }  of  a  PRS  with 
period  P  is  also  periodic  with  period  P 


Rp(k) 


(k)modP  ■  0 

o.w. 


(5.3) 


Consider  a  signal  s(n)  made  up  of  two  independent 
PRS's,  Pi(n)  and  P2(n)  with  periods  Pi  and  P2  respectively, 
that  is  s(n)Bpi  (n)'^p2(n) .  Assume  that  Pi  and  P2  are 
relatively  prime  and  that  our  observation  consists  of 
L*Pi.p2  samples  and  consider  the  detection  problem  in  which 
s(n)  is  received  in  additive  noise.  Now  the  P1.P2  samples 
in  the  observation  are  really  generated  by  P1+P2  independent 
samples.  This  situation  is  analogous  to  that  of  the  M.N 
peaks  in  the  ambiguity  diagram  being  generated  from  M-^N 
independent  parameter  sets.  In  Sect.  5.3  we  will  study  the 
detection  of  s(n)  (as  described  above)  in  noise  based  on  two 
design  hypotheses,  one  based  on  Pi. Pa  independent  samples 
and  the  other  based  on  Pi-^Pa  independent  samples.  Now  we 
study  the  detection  of  Periodic  Random  Sequences  in  noise. 


5.2  Detection  of  one  PRS  of  Known  Period  in  Gaussian  Noise 
Let  the  observation  be  a  vector  ^  of  length  L, 

rh 

statistics  of  jr  under  Ho  and  the 


61 


Hi  hypotheses  are 

Ho  :  ■  n(i)  ,  n(i)  “  N(0,1)  i.i.d.  (5.4a) 

Hi  :  y^^  ■  Pi(i)  ^  n(i)  ,  £i  and  n  independent  (5.4b) 

where  Pi(i)  is  a  PRS  with  period  Pi  and  pi (i )~N(0,Ai ) .  We 
assume  Pi  divides  L,  i.e.  Li.PibL,  where  Li  is  the  number 
of  periods  of  Pi(i)  in  the  observation.  Let  R  be  the 
autocorrelation  matrix  of  the  observation  under  the 

signal+noise  hypothesis,  R»E(^.y  |Hi).  For  now  assume  that 

R'^  exists,  so  we  may  write  down  the  probability  density 
function  of  the  observation  under  the  two  hypotheses 

f(y|H.,)  »  (5^)^'^^.exp(-y'^y/2)  (5.5a) 

(5.5b) 

where  |R|  is  the  determinant  of  R.  From  eq.  5.5  it  follows 
that  the  likelihood  ratio  &(y)  and  the  loy-likelihhod  ratio 
z(y)  of  the  observation  are  given  by 

a(xlB)  «  (j^)^/^.exp(^'^(I-R"’)y/2)  (5.6a) 

z<y|B)  -  ^z'^(I-R‘^)I  -  5ln(lR|)  (5.6b) 

where  0  {Ai,Pi}.  The  likelihood  and  log*- likelihood  ratios 
are  conditional  to  the  period  and  the  signal  power. 

Before  proceeding  further  we  need  to  take  a  closer  look 
at  the  autocorrelation  matrix  R 

R  -  E(y.y'^lHi) 


(5.7) 


Where  is  the  noise  autocorrelation  matrix  and  R_  is  the 
n  Pi 

autocorrelation  matrix  of  PRS  pi,  the  above  decomposition 

follows  from  the  independence  of  noise  and  PRS  Pi.  Now 

Rj^  -  Ij.  ,  the  L\L  identity  (5.8a) 

Rp^  ■  Ai{6( I i-j |modPi)}  i, j»0, 1 , . . . rL-1  (5.8b) 


Of  course  R  is  Hermitian  symmetric,  and,  as  such  is  similar 
to  a  diagonal  matrix  [31  pp.  201-202].  Due  to  PRS  Pi,  R 
also  has  the  interesting  property  that  each  row  of  R  is  a 
right  circular  shift  by  1  of  the  row  immediately  above  it. 
Matrices  with  this  property  are  called  circulant  or  cyclic 
matrices  [32,33  pp.  133-139]. 

5.2.1  Properties  of  Circulant  Matrices 

Since  the  theory  developed  in  this  chapter  relies 
heavily  on  the  properties  of  circulant  matrices,  some  of  the 
useful  properties  are  summarized  here.  Let  be  a  nxn 
circulant  matrix,  then 


C 


n 


Co 


Cl 

Co 


Ca.  ..c  1 

ci...c„.2 

Co« • •c^.3 


•  4 


(5.9) 


Circulant  matrices  are  a  special  case  of  Toeplitz 
matrices.  The  matrix  may  easily  be  diagonalized.  Let  1(1^  be 
the  rth  eigenvalue  and  u^  the  rth  eigenvector,  then  and 
U|.  are  the  solutions  of 
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or  equivalently  of  the  following  system  of  difference 
equations 


m-1  n-1 

^'n-iii+k“rk  *  ^'k-m^rk 
k»5  k«m 


f  "  0,l,...,n-1(5.10b) 
^  r  -  0J,...,n-1 


where  4*^  is  the  rth  eigenvalue  and  is  the  kth  element  of 
the  rth  eigenvector.  It  is  easily  verified  [32,33  pp.  133- 
139]  that 


n-T 

ip  •  ^^Cjj.exp(“j2irrk/n)  ,  r  »  0,1..., n-1 


(5.11a) 


-r  “  (^)^^^.(1»exp<’j2irr/n),...,exp(-j2iT(n-1)r/n))'^  (5.11b) 
r  «  0, 1 , . . ♦ ,n-1 

The  sequence  of  eigenvalues  io*ii r •  •  •  is  the  Discrete 
Fourier  Transform  of  the  first  row  of  the  matrix  and  the 
eigenvectors  are  independent  of  the  elements  of  Now 
define  the  matrices 


(5.12a) 


(5.12b) 


The  matrix  is  a  unitary  matrix,  so  U*  «  conjtuj]  •  u”^. 
And  C|^  is  unitarily  similar  to  the  diagonal  matrix  that 
is  is  a  normal  matrix  [31  pp.  201-202]. 


n  n  n  n 


(5.13a) 


n  n  n  n 


(5.13b) 
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C~^  will  exist  if  and  only  if  all  the  eigenvalues  are  non¬ 
zero. 

5.2.2  Eigenvalue  Eigenvector  Decomposition  of  R 

Let  R(n)  be  the  first  row  of  the  matrix  R  ,then 
R(n)  »  6(n)  +  Ai6( {n)modPi)  ,  n  »  0,1,...,L“1  (5.14) 


Li-1 

-  6(n)  +  Ai^S(n-lPi)  ,  Li  » 

The  sequence  of  eigenvalues  i(>(m)  of  R  is  given  by 


t)i(m) 


L-1 

y"'  R(n). 
n»0 


exp(-j2irnm/L) 


(5.15) 


L-1 

6(n)  .exp(-j2irnm/L)  + 

L-1 

Ai)  }  4(n-lPi).exp(-j2wna/L) 
n^  i-O 

^  I  1  ^AjLj  f  ro*  0 f L I y 2L ) ^ ^ (P )  —  1 ) L } 
t  1  y  O.W 

Let  •  diag(i|»{0),i(1)y...yi(i(L-1)),  be  the  diagonal 

matrix  of  the  eigenvalues  of  R  and  the  unitary  matrix  of 
eigenvectors  of  R.  Then 


R 


vX 


(5.16) 


and  since  all  eigenvalues  are  non  zero 


5,2,3  The  Opt imum  Detector 

We  had  derived  the  log-likelihhod  ratio  z(£|3)  in 
eg,  5,6b,  Substituting  for  jR|  and  in  eg,  5,6b  and 

using  the  fact  that  we  obtain 

z(2|P)  -  “  llnl'i'Ll  (5.18a) 

where  •  Det  •  Det  R  -  (1  +  (5.18b) 

Denote  the  diagonal  matrix  by  6^^  and  let  e(m) 

denote  the  seguence  of  the  diagonal  entries  of  Q^,  8(m)  is 

then  simply  given  by 

9(»)  •  '  *  *Tm) 

,  I  AiLi/(  l+AiLj )  »  ro  ■  0(Li  ^ » r  (Pi^l  )1^i  (5,19) 

t  0  o.w. 

For  brevity  define  Ci  ^  AiLi/( HAiLj ) 

So  now  z()r|6)  is  given  by 

2<ZlB)  -  <^*20) 

Now  there  are  two  options  available  to  us  1)  Calculate 
the  vectors  Uj^  and  first  or  2)  Calculate  the  matrix 

consider  both  options,  the  first 
option  gives  us  better  insight  and  the  second  leads  to  a 


practical  detector  structure 


5.2«3. 1  Optimum  Detector  in  the  Frequency  Domain 
Define  the  row  vector  Y  as  follows 

^  (5.21a) 

it  follows  that 


(e)’^V  ■  tjjz  ‘S-S'b* 

The  row  vector  Y  is  simply  the  DFT  of  2.  with  y  real,  Y* 
is  the  conjugate  transpose  of  Y.  So  zC^lp)  is  equivalent  to 
the  statisitc  z(Y|B)  given  by 

(5.22) 

If  we  assume  that  the  period  Pi  and  the  power  level  Ai 
are  known  an  equivalent  detection  statistic  is  ^(Y) 


♦  (Y|0)  ^  *  (5.23) 

Since  is  a  diagonal  matrix,  we  may  simplify  the  above 
equation  as  follows,  let  Y(m)  be  the  mth  element  of  Y  then 


♦(IlB) 


1 


L-1 


0(m) I Y(m) I  * 


(5.24a) 


m 


I  Y(m)Y(m) I  * 


(5.24b) 


where  y{m)  ^  (^0(m))^' 


(5.24c) 


Detection  based  on  4(Xl5)  is  similar  to  an  energy 
detector,  we  will  call  it  a  "  frequency  selective  energy 
detector  "  or  FSED.  In  fact  yim)  defines  a  digital  comb 
filter  [34  p,  241],  A  block  diagram  of  the  above  detector 
is  given  in  fig.  5.1.  We  now  examine  the  sequence  Y(m) 
under  both  the  noise  alone  and  the  signal  +  noise 
hypotheses. 


Ho  :  y(i)  “  n(i)  ,  i  *  0,1,,..,L“1 

&  y(m}  a  N{m)  ,  u\  a  0,1,,.,,L-1 

Hi  :  y(i)  a  pj(i)  +  n(i)  ,  i  a  0,1, ...L-1 

&  Y(m)  a  Pi{m)  +  N(m)  ,  m  a  o,1,...,L-1 


(5.26a) 


(5.26b) 


Where  N(m)  is  the  L  point  OFT  of  the  observed  noise  sequence 

and  Pi(m)  is  the  L  point  DFT  of  the  PRS  pi(i).  Define  pi(i) 
to  be  one  period  of  PRS  pi(i). 


Pi(i) 


Pi(i)  ,  i  a  0,1,...,Pi-1 
0  ,  o.w. 


(5.27) 


then 


Pi  (m) 


L-1 

Fp.(i 

!■() 


)  .exp(-j2iiim/L) 


(5.28) 


Li-1  (n+1)P,-1 


a  Pi(i-nPi).exp(-j2Trim/L) 

nafl  i»nPi 


and  usina  the  substitution  s  •  l-nPt.  we  oet 


diagram  of  the  Frequency  Selective  Energy  Detector 
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Li-1  Pri 

^  ^  Pi(s?  .exp(~j2tr(s+nPi)m/L) 

na6  S’s'O 


Li-1 

now  ^  exp{-j2irnm/Li) 
n»0 


Lj  ,  ni  ■  0 f L j  f  •  •  • ,  (P} *“  1  )L ) 

0  ,  o.w. 


So 


Pi-1 

I  pi(s).exp(-j2Trsm/L)  ,  in«0,L j , . .  * ,  {Pj^ t  )Li 

Pi(m)  -  i  (5.29) 

0  f  o.w. 

The  above  shows  that  under  H}.,  Y(m)  has  signal 
components  only  for  , . , ♦ , (P|~1 )Lj .  This  is 

intuitively  pleasing  because  0{sO  Is  non-zero  for  precisely 
the  above  m's.  So  the  f^e^^uensy  selective  energy  detector 
rejects  the  *out  of  band  noise." 

5. 2. 3. 2  Opt imum  Detector  in  the  Time  Domain 

Now  we  first  calculate  the  matrix  will 

lead  to  a  time  oomain  solution  for  the  detector.  Let 

Q  ^  t'nen  because  6]^  is  diagonal  Q  is  circulant.  So 

the  matrix  Q  may  be  determined  in  terms  of  its  first  row 
Q(n).  in  fact,  the  first  row  is  just  the  inverse  DPT  of  the 
sequence  e<m}. 


v>'>  /•IT* V.**  iy  ,.** t**. .*■*  .**,-■*, 


L-1 


Q(n)  a  ^  ^  ©(a) .exp(  j2irnin/L) 


L“1 


Pi-1 


(ci/L). 


mcO  r»0 


6(m-rLi)  .eKp(  j27rnm/L) 


(5.30) 


Pi-1 

s  (cj/D.)  exp( j2irnr/Iii) 
r*0 

s  ;  c  1 1*  1  /L  ,  n  “  C»Pi,...,(Li”'1)Pi 
:  0  ,  o.w. 

So  2(^1 0)  is  given  by 


(5.31a) 


Pi-ILi-1 


(5.31b) 


**  *■’  4*aiW|P 

3^  KaO 


If  we  assume  that  the  signal  power  and  period  are  known 
then  an  equivalent  detection  statistic  is  ^(£)  giver  by 


♦(2le) 


fCjPi 

"3r 


Pi-ILi-1 

.^{^y(j^kP»))’ 

3t*iy  k«o 


(5.32) 


We  call  the  detector  that  uses  ^(£l5)  as  the  detection 
statistic  the  circulating  average  energy  detector  (CAED)  or 
the  periodic  averaging  energy  detector.  Pig,  5.2a  is  a 
block  diagram  of  the  generic  signal  processor  that  produces 
^(yU)  with  the  parameters  P^,  L|  and  A(.  The  ''analog  shift 
register”  is  initially  set  to  zero,  once  the  desired  number 
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o£  periods  (Li)  have  been  measured  the  switch  feeds  the 
output  of  the  "analog  shift  register"  into  the  squarer.  In 
future  we  will  call  this  signal  processor  the  CAEDSP  with 
the  appropriate  parameter  values  denoted  in  brackets. 
Fig.  5.2b  is  a  block  diagram  of  the  above  CAED. 

5.2.4  Performance  Evaluation  of  the  Optimum  Detector 

First  of  all  note  that  4>(2|g)  and  4(Ylp)  are  the  same 
random  variable,  hence  decisions  based  on  ()>(£|3)  are 
equivalent  to  decisions  based  on  $(Y|g)  and  they  result  in 
identical  performance.  In  this  section  we  derive  the 
probability  density  function  for  (j)(x|E)  and  z(^|g)  and 
compare  the  performance  of  the  optimum  detector  with  that  of 
the  total  energy  detector. 


5.2.4. 1  Derivation  of  the  p.d.f .  of  (i)(y)  and  z(y  I B) 

Lj-I 

Let  a(j)  “  y(j+kPi),  then  a(j)  is 

normal.  The  statistics  of  a(j)  are  derived  below 


E[a(j)lH,^]  »  0  ,  k  -  0,1 


(5.33a) 


E[a  (j)|Ho] 


Li“1Li-1 

y” E[y(j+kPi)y(j+sPi)|Ho]  (5.33b) 


(^1- 


6(  (k-s)Pi ) 


CxPiLi/2L  ■  c i/2  ^  Oo 


Fig,  5.2b  Block  diagram  of  the  Circulating 
Average  Energy  Detector. 


Li-lLj-l 

ECa*(j)|Hi]  -  (1^).^  ^E[y(j+kPi)y(j+sPi)  |Hi]  (5.33c) 


“  ^E[n( j+kPi)n( j+sPi)  +  p!(j)] 


2L 


+  AjLi  ) 


(Ci/2).(1  +  AjLi)  -  AiLi/2  ^  0? 


So  we  have 


£(a(j)|H|j)  =  (5j)’''^.(5-).exp(-  a^)  ,  k=0,t  (5.34) 

k  20|, 


Now  let  Y^j)  “  a  then  it  follows  that 


f(Y(j)|H,,) 


( - 5 ,  t(3)>0 


SttyC  j)a)j 


20i 


(5.35) 


,  o.w, 


which  is  a  gamma  density  function  with  one  degree  of 
freedom.  Now  the  yU)’s  are  independent  and  identically 

Pi-1 

distributed  under  both  hypotheses.  Since  (^(^IS)  »  ^  Y^dK 

<t)(y|6)  also  has  a  gamma  density  but  with  Pi  degrees  of 
freedom. 


f  (4.|H,^) 


,(Pi/2)-1 

(2oi[)^‘^^.r(Pi/2)  2o 


'  .exp(-  ,  4»  >  0 


(5.36) 


,  o.w, 


Where  r(.)  is  the  gamma  function  [35  pp.  255-263].  Now 
z(y|3)  -  ^♦(y|B)  -  •j.lnjfj^l ,  substituting  for  gives 
z(y|B)  -  ^4>(z|p)  -  (Pi/2).ln(l+AiLi).  If  we  let 
bi6(Pj/2) .ln( 1+AiLi )  then  it  follows  that 


f(zlHjj) 


(20j^) 


r(Pi/2) 


2oi 


(5.37) 


,  o.w, 


If  u  is  a  threshold  level  for  z,  the  probabilities  of 
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detection  and  false  alarm  are  given  by 


Pj3  =  ;*f(z|Hi)dz 

(z+bi)<P^/2)-1  z+b, 

^  /•«  _ .exp( - r)dz 

(2ot)‘’‘^^.r(p,/2) 


using  the  substitution  w  «  z+b^,  we  obtain 


Pd  "  ^ 


CD 

u+b 


(Pi/2)-1 

.exp( - ^V)dw 

1  “^^772 - 


(2oi)*^‘'‘‘.r(Pi/2) 


2o  1 


Pi  u-»-bi 


r(3 


)/r(p,/2) 


2o  1 


1  -  Y(tL.  “*';‘l/r(p./2) 
^  2^ 


similarly 

PpA  •  r(^.  ^)/r(p./2> 

‘  2a, 

-  '  -  y(^.  :!!^)/r(p,/2) 

^  2a, 

via,  a)  and  Y(.r>)  are  the  incomplete  gamma 
defined  as  [35  pp.  255-263] 

r(a,x)  ■  /“  t'^'^expC-tidt 


(5.38a) 

(5.38b) 

(5.38c) 

(5.38d) 

(5.39a) 

(5.39b) 

functions 

(5.40a) 


yia,yi) 


/q  t®^”^  .exp(-t)dt 


(5.40b) 
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5, 2, 4, 2  R,O.C.  Curves  and  Comparison  with  Energy  Detector 

The  energy  detector  is  a  sub-optimum  detector,  it  would 
be  optimum  if  the  signal  samples  were  independent.  That  is, 
it  ignores  the  signal  periodicity.  So  the  energy  detector 
is  based  on  the  following  observation  statistics 

Ho  =  y^  ~  N(0,1)  i.i.d.  (5.41a) 

HI  »  y^  **  N(0,l+Ai)  i.i.d.  (5.41b) 

and 

f(y;|Ho)  “  (5^)^'^^.exp(-  ^^/2)  (5.42a) 

£(Z|HU  - 


The  likelihood  ratio  and  the  log-likelihood 

ratio  Zg^(l)  are  given  by 


(5.43a) 


A 

"  (L/2).ln(1+Ai)  (5.43b) 


is 


the  total  energy  in  the  observation. 


Assuming 


signal  power  is  known  an  equivalent  detection  statistic  is 
•n(^)  given  by 


Tl(l) 


Ai 

ITT+aTT 


L-1 


(i) 


(5.44) 
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A  block  diagram  for  the  energy  detector  is  given  in 
fig.  5.3.  Under  Ho»  and  have  gamma  density 

A  2 

functions.  If  we  let  b2=(L/2) .ln( 1+Ai ) ,  and  with  Oo=l  we 
have 


f (nlHo) 


.exp(-  -^)  ,  i\  >  0 


o.w. 


(5.45a) 


f(Zgdl^o)  “ 


(2oo}^/^r(L/2) 


0  ,  o.w. 


T 

2oo 


(5.45b) 


If  eq.  5.41b  were  true,  then  under  HJ,  t\(^)  and  would 

also  have  density  functions  similar  to  eqs.  5.45a  and  5.45b 

2  2 

respectively;  with  Oo  replaced  by  oi-l+Aj.  If  we  assume 
that  eq,  5.4b  is  true,  then  the  density  function  of 
Zgd^Z)  cannot  easily  be  found  because  the  terms  in  the 
summations  of  eq.  5.43b  and  5.44  are  not  independent. 
However,  it  is  still  possible  to  compute  some  relevant 
statistics  of  'n(z)  and  z^^(^)  under  and  the  assumption 
that  eq.  5.4b  is  true. 

In  fig.  5.4a  the  R.O.C.  curves  for  the  CAED  are  given 
for  several  signal  power  levels.  Noise  is  assumed  N(0,1). 
We  assume  L-1000,  the  signal  periods  range  from  21  to  30 
inclusive.  The  number  of  periods  measured  in  an  observation 
is  the  number  of  integer  periods  in  a  sample  size  of  1000. 
We  have  plotted  the  average  probability  of  detection  vs. 
the  probability  of  false  alarm.  The  average  detection 
probability  is  obtained  by  finding  the  detection  probability 
for  each  signal  period  for  a  given  false  alarm  probability 
and  then  calculating  the  average  of  the  detection 
probabilities.  For  comparison  purposes  in  fig,  5.4b  the 
R.O.C.  curves  for  the  ED  are  also  given.  The  R.O.C. 


.55  .95 


for  signal 
16,  .032. 


Pig*  5.4b  R.O.C.  curves  for  the  ED  with  the  same  signal 
and  noise  statistics  as  in  fig.  5.4a. 
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curves  are  obtained  using  Monte  Carlo  simulation  methods. 
The  details  of  the  simulation  are  given  in  Appendix  A. 

Analytically  we  can  compare  the  performance  of  the  two 
detectors  using  the  concept  of  normal  detectability.  As 
given  in  Chapter  II  the  definition  of  normal  detectability 
is  valid  for  simple  gaussian  hypotheses.  Here  we  assume  the 
gamma  densities  of  eqs.  5.36  and  5.44  are  approximately 
gaussian.  This  will  allow  us  to  obtain  approximate 
performance  figures.  The  expression  for  normal  detectabilty 
'd'  is  repeated  below 

d  ^  [  E(z|Hi)  -  B(z|Ho)  (5.46) 

vaFTzT 

As  defined  d  is  a  measure  of  the  output  SNR,  We  will 

look  at  both  d  and  d*  »  (d)^^^.  For  simple  gaussian 
hypotheses  var(z)  is  the  same  under  each  hypothesis.  This 
is  clearly  not  the  case  here.  However  if  we  assume  that 
signal  power  is  small  then  we  can  use  the  approximation 
var(z|H|)  a  var(z|Ho)  [16],  Because  ^  differs  from  z  by  a 
constant  we  need  only  calculate  E[4|H|^]  ,  k  «  0,1  and 

var(4),  similarly  with  the  energy  detector.  For  the 
normalized  gamma  density  function  f(x)  the  first  and  second 
moments  are  given  by 


f(x)  ■  |i^.exp(-x)  ,  0<x<» 

(5.47a) 

®  xP 

“  h  rXp7*eJ‘P<'“x)dx  -  p 

(5.47b) 

2  00  xP'*’^  2 

E(x  )  -  /q  i^.exp(-x)dx  -  p  ♦p 

(5.47c) 

Thu  first  and  second  moments  of  the  CAEO  detection 
^  are  now  easily  obtained. 

statistic 

e4|h,^)  »  /q  4..f(4>|Hj^)d<^ 


using  the  substitution  x  ■  A,  we  obtain 

2ot^ 


2 


E(*lH|j)  •  20|j./q  j|p-7jy.exp(-x)dx 


2 

'k* 


Oi.Pl 


E(^  IH,^)  ■  /J  ♦  f(^|H,j)d* 


using  the  same  substitution  as  above  we  obtain 


a. 2  ,<P./2)^1 


E(4>  iHjj)  •  (20|^)  ./q  — .exp(-x)dx 


-  <2o’)\(^  > 


J .  a 


and  var(^(H|j)  »  (2a^)  .Pi/2 

Prom  the  above  we  can  easily  write  the  statistics 
and  Hi 

E(^|Ho)  •  Oo*Pj 


(5.48a) 


(5.48b} 


(5.48c) 


(5.4Sd} 


under  Ho 

(5.49a) 


CiPi/2 
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«  o\,?i  .  (5.49b) 

=  CxPi.(l+AiLi)/2 

var(4»|Ho)  «  (2o5)\pi/2  (5.49c) 

a 

=Ci.Pi/2 

var(4>|Hi)  a  (2o!)^.Pi/2  (5.49d) 

■  Ci.Pid+AjLj)  V2 


can  now  obtain  an  expression  for 
above  in  eq.  5.46 


[c,Pt(UA,Li)  >  CiPi]^ 
- 


d 


by  substituting  the 


(5.50) 


•  Ai*L.Li/2 


Going  through  a  similar  set  of  calculations  for  the  energy 
detector,  we  find 


E(iilHj)  -  AiL/2 
E(iilHa)  -  AiL.  t 

ITTOTT 

var(n|Ho)  •  A?L.  1 

^  (UA|) 


2 


(S.Sla) 

(5.51b) 

(5.51c) 


\  •  aJ.L/2 


(5.52) 


For  fixed  performance  ,  i.e.  d^  >  d.  and  with 

T]  9 
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»  input  SNR  for  the  ED  and  *  input  SNR  for  the  CAED  we 
have 


A^.L  =  A^.L.Li  (5.53) 

^  1  *  I 


or 


itt  ■  (5.54) 

% 

Comparing  d  for  the  CAED  and  the  ED  for  fixed  input  SNR  we 
obtain 

S  .  L,  <5-55) 

n 

The  CARD  outperforms  the  ED  by  lOlogLt  dB*;  for  our  case 
ranges  from  33  for  Pt«30  to  47  for  Pi»21  so  the  CAED 
performance  is  approximately  IS. 2  dB  to  16.8  dB  better  than 
the  ED  performance.  For  fixed  input  SNR  the  CAED 

performance  goes  up  linearly  with  This  approximate 

result  is  in  reasonable  agreement  with  the  simulation 
results.  (Reminder  :  Lt  is  the  number  of  periods  ^looked 
at**  in  one  observation.) 

5.2.5  Detection  of  PRS* a  With  Unknown  Period 


A  related  problem  of  interest  and  in  some  sense 
analogous  to  the  classical  "Defection  of  one  of  M  Orthogonal 
signals’*  problem  is  the  problem  of  detecting  a  PRS  in  noise 
where  the  period  is  one  of  M  possible  periods.  The 
observation  under  the  two  hypotheses  Ho  and  Hi  is 


Ho  :  y(i)  *  nil)  ~  N(0,1)  i.i.d 


(5.56a) 
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Hi  :  y(i)  =p(i)  +n(i)  ,2&n  independent  (5.56b) 

where  p(i)  is  a  PRS  as  described  in  Sect.  5.1  with  period  P 
and  where  P  is  an  element  of  the  set  S2«{Pi  ,P2 » .  • .  f .  We 


assume  that  P=P|j  with  probability  Pr(Pj^).  We  let  the  number 
of  samples  in  the  observation  equal  L.  We  define  Lj^  as 
follows 

^  I  L/Pj^  I  (5.57) 

where  l_x_(  »  largest  integer  no  larger  than  x.  So  that  Lj^ 
is  the  number  of  integer  periods  of  PRS  with  period  Pj^  in  an 
observation  L  samples  Icng.  We  assume  for  simplicity  that 
all  sequences  have  the  same  power  level  A,  We  define  the 
log-likelihood  ratio  given  PRS  with  period  Pj^  as  follows 

A 

where  Cj^ 

In  all  cases  the  optimum  detection  statistic  is  the  “average 
likelihood  ratio”  Following  from  the  Bayes-Birdsall  Theorem, 
let  g  be  a  vector  of  signal  parameters  with  known  joint 
distribution  and  density  fi^),  let  y  be  the  observation  and 
f(j|ly)  be  the  a-posteriori  density  of  then 

t(2)  -  tlylB).^‘SI  (5.S9a) 

nuF 

it  follows  that 


-  ^k.ln(HAL.)  (5.58) 

-Hjr  K 


£(y)nfily)  *  A(ylj|)f(fi) 


(5.59b) 
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Integrating  both  sides  with  respect  to  £ 

!Hz)  »  /  £(2li)f(i)di  (5.60) 

For  one  of  M  periods  we  obtain 


M 


iiz)  -  )  £(z|P)t)*Pr(P|i) 


(5.61) 


If  we  assume  equally  likely  to  be  any  member  of  $2  we 
obtain 


(5.62a) 


M 

«  ^  ^  exp(z(jr  jP|^) )  (5.62b) 

fcrr 

M 

z(^)  »  ln(^^exp(z(y|Pj^) ) )  -  InM  (5.62c) 


A  block  diagram  of  the  optimum  detector  is  given  in 
fig.  5.5.  The  R.O.C.  curves  are  given  in  fig.  5.6.,  for  (I 
*  (21 ,22, . . . ,30}  and  L«1000.  Noise  as  usual  is  assumed 
N(0,1)  i.i.d. 

Next  we  compare  the  performance  of  the  above  detector 
with  that  of  the  estimator-detector.  The  estimator-detector 
calculates  z()!^|Pj^)  for  all  P^  c  R  and  bases  its  decision  on 
max  ^(xlPi)*  I'he  block  diagram  for  the  estimator-detector 

Pj€Q  ^ 

is  given  in  fig.  5.7  and  the  R.O.C.  curves  in  fig.  5.8. 
For  comparison  purposes  the  R.O.C.  curves  for  the  energy 
detector  are  also  given  in  fig  5.9.  From  the  R.O.C.  curves 


-Bln(l+A.L.) 
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Block  diagraa  of  the  optimum  detector  when  the  signal  period 
is  equally  likely  to  be  one  of  M  periods. 
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K 


Pig.  S.6  Averaged  R.O.C.  curves  for  the  optimum  detector 
of  fig.  5.5  vith  M«10  and  signal  powers  of  *002. 
.004.  .008.  .016.  .032  and  .064. 


^  Ca£OSP(f>k,Ai,U 
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Pig.  5.7  Block  diagram  of  the  estimator-detector  when  the  signal 
period  is  equally  likely  to  be  one  of  M  periods. 


Fig,  5,8  Averaged  R.O.C,  curves  for  the  estimator-detector 
of  fig,  5.7  with  M«10  and  signal  powers  of  .002, 
.004,  .008,  .016,  .032  and  .064. 


Fig.  S.9  R.O.C.  curves  for  the  energy  detector  with  signal 
powers  of  .002,  .004,  .006,  .016,  .032  and  .064. 


for  the  optimum  detector  and  the  estimator-detector  we  see 
that  the  optimum  (average  likelihood  ratio)  detector  and  the 
estimator  detector  have  essentially  the  same  performance. 
Monte  Carlo  methods  were  used  for  generating  the  R.O.C, 
curves. 

Note  that  the  P.  for  which  z(y|P,*)  is  maximum,  is  the 
maximum  likelihood  estimate  of  the  period  of  the  PRS.  Plots 
of  Pr(  Error  in  Estimation  of  period)  vs.  signal  power  are 
given  in  fig.  5.10.  We  now  look  at  zCilPj^)  in  more  detail. 


Pi-1  L.-1 


•  i'-rSb-^'S 


-  ^i.ln( 1  +  AL, ) 
”7  i 


(5.63a) 


Pi-I^l 


L-l 


1. 

1  r 


-  *^i.ln(1  +  AL,)  (5.63b) 

-j  i 

\  .  * 

but  A  y(j’''<nP^)  is  the  m.m.s.e*  estimate  of  p^(j) 

(IS  pp.  56-59],  the  jth  sequence  sample  under  the  assumption 
that  the  sequence  period  is  P^;  and  where  0  <  j  <  P^-1.  We 

A 

define  p^(j)  as 


L.-1 


Pit])  ^  A  .1  .y~~ y(j^mP.) 

m»0 


(5.64) 


Fig.  5.10  Average  Probability  of  Error  in  estimating 
signal  period  vs.  S/H  ratio. 
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Then  it  follows  that 


Y~ y^3+kPi)*Pi(j) 
^0  K»C 


^i.lnd 


+  AL. ) (5.65a) 


or 


z(£lP^)  »  “  ^i.lnd  +  AL^) 

"2 


(5.65b) 


where  2i  estimated  PRS  waveform  under  the  assumption 
that  the  PRS  period  is  P^.  The  detection  statistic  may  now 
be  written  as 


max  z(y |P. ) 
“i 


max 

"i 


in  ^ 

{  ‘Ej 


Ind  A.L^  } 


(5.66) 


Because  of  the  simplicity  of  the  structure  of  the 
estimator-detector  (  no  exponential  non-linearity)  and  no 
loss  of  performance  compared  to  the  optimum  detector,  the 
estimator-detector  is  the  preferred  realization  of  the 
optimum  detector  in  this  case.  The  advantage  of  the 
estimator-detector  over  the  energy  detector  are  1)  Better 
performance  2)  The  ability  to  simultaneously  estimate  the 
signal  period. 

The  problem  of  estimating  signal  period  (or  frequency) 
is  an  interesting  problem  in  its  own  right  [36,37].  The 
above  development  focused  on  a  periodic  random  sequence.  We 
now  outline  the  development  for  deterministic  but  unknown 
periodic  sequences.  This  development  differs  from  the 
previous  in  that  it  is  non-Bayesian;  the  signals  are 
•unknown",  not  random  vectors. 

Assume  the  observation  consists  of  L  samples.  The 
statistics  of  the  observation  are 
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Ho  :  y(i)  ■  n(i)  **  N(0,1)  i.i.d. 


(5.67a) 


Hi  :  y(i)  »  n(i)  + 


(5.67b) 


where  is  a  periodic  sequence  of  period  ^  * 

We  assume  that  the  period  P|^  is  unknown. 

Let  Lj^  »  snd  define  Pjj(i)  as  follows 


m0 

Pk 


Pjj  ( 1 )  y  i  "0 » 1  f . . . » Pjj"“  1 

0  t  o.w. 


(5.68) 


If  the  sequence  p,.  is  known  then  the  optimum  estimate 
of  the  period  follows  from  forming  the  M  log-likelihood 
ratios  z(y|P|^)*  k«1,2,...,M  and  choosing  the  largest. 


f(ylp^)  -  (yj)  ‘^.exp(-  5(l*fi|j)'^.(rEitM 

(5.69a) 

f(y|Ho)  - 

(5.69b) 

and 

2(xlPjt)  - 

(5.70) 

However  as  the  sequences  are  not  known  we  follow 
custom  and  first  form  a  maximu’^  likelihood  estimate  of  the  M 
possible  sequences  based  on  the  observation,  and  then  use 
these  estimates  in  eq.  5.70  and  then  choose  the  maximum.  We 
can  rewrite  2(y|Pj^)  as  follows 

Vk“’ 

z(y!Pjj)  “  (y(i)  -  5P,j(i)).pj^(i) 

(5.71a) 
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(^(y(j+rPji) 

j»0  r»0 


(5.71b) 


Differentiating  eq.  5.71b  with  respect  to  each  P|j(j)  and 
setting  the  result  to  zero  we  obtain 


-  ^(y(j^rP,j)  -  p,^(j))  -  0  (5.72) 

3P|j(j)  ^ 


for  3  * 


or 


(5.73) 


Pl^(j)  is  the  oiaximum-liRelihood  estimate  of  P|iU)* 

A 

A  «» 

sequence  Pj^  extends  p|^  periodically 


Pi,(i) 


V. 

r-“  - 

L  '’k' 


i“rPj^)  0  i *0 » 1  * . . . « Ii“  1 


.  o.v. 


(5.74) 


Based  on  the  estimated  sequences  we  form  the  M  statistics 

A 

z(^|P|j)  defined  as 

z(xlPjj)  ^  <2  “  1*£r^^‘Er  (5.75) 


The  estimated  period  is  then  the  one  that  corresponds  to 
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max  z(£|pj^).  The  block  diagram  of  the  estimator-estimator 
Pk 

is  given  in  fig,  5.11. 

5.3  Detection  of  two  or  more  PRS*s  in  Noise 

In  this  section  we  generalize  the  results  obtained  in 
Sect.  5.2.  We  first  examine  the  problem  of  detecting  2  or 
more  PRS's  with  known  periods.  Next  we  examine  the  problem 
of  the  detection  of  k  of  M  PRS’s  in  noise.  Finally  we 
develop  the  equations  for  the  most  general  problem  that  fits 
into  the  framework  established  here,  in  this  case  we  allow 
signal  and  noise  to  be  complex  with  the  noise  not 
necessarily  white  and  successive  signal  samples  not 
necessarily  independent. 

5.3.1.  Detection  of^  ^  or  more  PRS* s  of  Known  Periods. 

In  this  case  the  signal  consists  of  the  sum  of  N  PRS's 
(N>2)  of  known  periods  Pt ,P2,. . The  number  of  samples 
in  the  observation  is  chosen  so  that  L»L,C.M,  {Pi  ,P2, . . .  ,Pj|} 
or  some  integer  multiple  of  the  L.C.M.  This  insures  that 
each  PRS  has  an  integer  number  of  periods  in  the 
observation.  We  define  Lj^  as  the  number  of  periods  in  the 
observation  of  PRS  with  period  Pj^  i.e.  ■  ^/^k 

The  observation  under  the  two  hypotheses  Hq  and  Hi  is 

Ho  :  y(i)  -  n(i)  "  N(0,1)  i.i.d.  ,  i-0, 1 , . . . ,L-1  (5.76a) 


N 

:  y(i)  ■  n(i)  +  ,  i«0, 1 , , . .  ,L-1 


(5.76b) 


Where  P)j(i)  has  period  Pj^  and  p|^( i)~N(0,Aj^) .  Noise  is 
assumed  independent  of  the  PRS’s  and  the  PRS’s  are  assumed 

independent  of  each  other.  We  let  R>E(y.^  |Hi)  and  write 


.•  1,  ■ .  'x  ’ji  • 
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the  density  function  of  the  observation  under  the  hypotheses 


Ho  end  Hi  as  follows 

(5.77a) 

fs’’®"'!)  (5.77b) 

The  likelihood  and  the  log-likelihood  ratios  follow 

il(y)  ■  (-j^)^'^^.exp(^.2'^(I-H’by)  (5.78a) 

2(y)  »  ^.y'^(I“R“bx  -  ^.ln|R|  {5.78b) 


As  before  the  key  to  the  simplification  of  the 
equations  for  illy;)  and  lies  in  the  eigenvalue- 
eigenvector  decomposition  of  R.  First  we  write  R  as 

R  •  E(y.y^|Hi)  (5.79) 


Where  R^^  is  the  noise  autocorrelation  and  is  the  LxL 
identity  matrix,  R^  is  the  autocorrelation  matrix  of  PRS  p|^ 
and  is  circulant.  The  above  decomposition  follows  from  the 
independence  of  the  noise  and  the  PRS's.  It  follows  that  R 
is  also  circulant.  In  fact 


R  «  I, 


N 

*  ^  A|j{6(|i-j|modPjj)) 

fcr 


.,L-1 


(5.80) 


The  first  row  of  R,  R(n)  is  given  by 
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N 


R(n)  »  6(n) 


AkY” 5(n-lPu)  ,  n=0,1,...,L-1 


t|>(in),  the  eigenvalue  sequence  of  R  is  given  by 


(5.81) 


L-1 

i|)(in)  «  )  R(n) .exp(-j2irnin/L) 
n»0 


(5.82a) 


N 


A|^L|^  I  in*0 


k«1 ,2, , . . ,N 
1  ,  o.w. 


(5.82b) 


Leu  ^  «  diag(<i(0)  ,tli(  1 ) , . . .  ,^(L-1  / )  be  the  diagonal  matrix  of 
eigenvalues  of  R  and  be  the  unitary  matrix  of 
eigenvectors  of  R,  then 


B  .  Uj_rj* 


(5.83) 


and  since  all  eigenvalues  are  non*-zero 


R'^  • 


(5.64) 


also 


|R|  •*  |y! 


(5.a5a) 


N  N 

•  (1  ♦  )  ^k^k^’  ” 
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N 


(1 


N 

n  (1  +  A.L. ) 
Jt-1 


N  ] 

..  n  (1  +  A.L. ) 
k»1  ^ 


(5.85b) 


Substituting  for  r"^  and  |R|  from  eq.  5.84  and  eq.  5.85  in 
eq.  5.78b,  we  obtain 


Z(l)  - 


N 


N 

(1  +  Y~  Aj^L,^) 


“  *10 


•& 


N 

n  (1  +  A.L.) 
k-1  ^ 


(5.86) 


We  define  0  ^  0  is  a  diagonal  matrix  and  the 

elements  along  the  diagonal  are  given  by  0(m)  •  1  -  1  ,  so 

7tm) 


0(m) 


£‘-‘‘ 


,  m*0 


N 


'•S 


AuL, 


^k^k 


*k"k  ,  m-L|^,2Lj^,...,(Pj^-1)Ljj 

^  ^k^k  k-1,2,...,N 


o.w. 


(5.87) 


Now  define  Co  and  C|^  as  follows 


102 


N 

■C.  H  ^ 


N 

_  \ 


k"lc 


N 


'  •  £*■■■■ 


Cu  ^  ^k^k 

We  may  now  rewrite  e(m)  as  follows 


(5.88) 


(5.89) 


e(m)  «  )  9jj(m) 


■& 


where 


0j^(m) 


-Co6(m)  ,  k-0 


Cj^.^  ~6(m~rLj^)  t  k*l|2f.»»N 


(5.90) 


(5.91) 


Now  let  Q^U^OU^  and  decompose  6  as  done  Cor  d(m)  so  that 


N 

e  A  y  9|j.  The  diagonal  elements  of  are  given  by  Q^ita), 

N 

Define  Q,^  ^  Vk^L  ®  "  /  ^k*  definition  the 

k"Ty 

matrices  Qj^  and  hence  Q  are  circulant,  so  we  may  obtain  Q  in 
terms  of  its  first  row.  Let  the  first  row  of  Q  he  Q(n), 
then 


Q(n/  ■  Qo(n)  +  Qi(n)  ...  ^  Qjj(n) 


(5.92) 
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and 

Qo(n)  » 

Q|^(n)  ■ 
It  follows 

z(x)  ■  • 


-  fjl  all  n 
L 


{5.93a) 


CuP 


k  k  /  n^O  ^ P|^ f  2Pj^ f  •  •  •  /  1 ) « P|j 


0  o.w. 


that 


(5.93b) 


N  N 


(1 


^In, 


N 

•£ 


"S 

n  (1  +  AfcLfc) 

k«t 


L”  1  N  ^k"*^  ^k”^ 

.{^y(n))*  +  ^°k^k.^  (^y(^rP,^)) 


Co 

IE 


N 


(5.94a) 


FO  bo 


(1 


^!k.ln(1  .  A|,L|,)  -  ^In 
kST  ^ 


N 

•£ 


n  (1  +  A^u) 

k-1  ^  ^ 


(5.94b) 


The  above  reduces  to 
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z(:^) 


N 

^z{z|P^) 


L-1 

“"(y y<n))* 


2l 


n»0 


N 


„  .  ^v>' 


N 

n  (1  +  A.Lfc) 

k-1 


(5.95) 


The  last  term  in  eq.  5.95  does  not  depend  on  the  reception, 
in  fact  it  is  a  constant  and  maybe  ignored  if  we  assume  we 
know  the  power  levels  Ay^'s.  The  second  term  in  eg.  5.95  is 
a  D.C.  correction  term.  It  is  at  D.C.  (zero  frequency) 
that  all  the  PRS’s  interact.  The  R.O.C.  curves  for 
detectors  based  on  z(y)  are  given  in  fig.  5.12,  for  N>2, 
Pib21,  P2*29  and  L«609.  Experimentally  we  have  shown  that 


N 

decisions  based  on  ^"~z(y|Py^)  are  essentially  equivalent  to 
decisions  based  on  z(y).  We  call  the  detector  that  bases 


its  decisions  on 


the  sum-detector . 


R.O.C. 


curves  for  the  sum-detector  are  given  in  fig.  5.13,  the 
signal  and  noise  statistics  used  in  generating  fig.  5.13  are 
the  same  as  those  used  for  generating  fig.  5.12.  For 
comparison  purposes  the  R.O.C.  curves  obtained  using  an 
energy  detector  are  given  in  fig.  5.14.  Since  the 
conditional  log-likelihood  ratios  z(y|P|^)  are  not 
independent,  it  is  not  easy  to  derive  the  density  function 
for  z(y).  However,  we  can  calculate  the  statistics  of  z(y) 
necessary  to  obtain  an  approximate  expression  for  the  normal 
detectability  d.  The  details  of  the  calculation  are  tedious 
and  are  given  in  Appendix  B.  In  calculating  the  statistics 
we  assume  that  K>2  and  that  both  the  PRS's  have  the  same 
power  A/2  and  that  Pi^Fa*  obtain  the  following 
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Fig.  5.12  R.O.C.  curves  for  the  optimum  detector  for  2 

PRS’s  of  known  period  in  noise  with  Pi*21,  Pa»29 

and  L«609.  Each  PRS  has  identical  power  levels 
of  .001,  .002,  .004,  .008,  .016  and  .032. 
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Pig.  5.14  R.O.C.  curves  for  the  energy  detector  with  the 
same  signal  periods  as  in  figs.  5.12  ft  S.13  and 
signal  power  levels  of  .002,  .004,  .006,  .016 
and  .032. 
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expression  for  d 

d  a  A*L(1  +  Li  +  L2)/2  (5.96) 

In  the  preceding  development  the  observation  size  was 
fixed  at  some  integer  multiple  of  the  L.C.M.  of 
{Pi , . . ,  ,Pjj} .  This  places  a  restriction  on  the  observer.  In 
the  following  this  restriction  is  relaxed.  The  number  of 
samples  in  the  observation,  L  ,  is  now  arbitrary.  It 
follows  that  there  are  Lj^  ■  1_L/Pj^__|  integer  periods,  of  PRS 
with  period  P|^  in  the  observation.  The  observation  under 
the  two  hypotheses  is 


:  y(i)  ■  n(i)  "  N(0,1)  i.i.d. 

(5.97a) 

N 

:  y(i)  *  n(i) 

(5.97b) 

The  Lxl«  autocorelation  matrix  is  no  longer 
circulant.  The  key  to  obtaining  the  equations  for  the 
optimum  detector  was  the  eigenvalue^eigenvector 

decomposition  of  the  circulant  matrix  R  and  R~\  For  the 
present  situation  where  the  number  of  samples  in  the 
observation  is  not  divisible  by  all  the  periods,  we  will 
call  R  the  expunged  circulant  matrix.  A  simple  well  defined 
procedure  for  the  eigenvalue-eigenvector  decomposition  of 
the  expunged  circulant  matrix  does  not  exist.  Also  a 
straight  forward  procedure  for  finding  the  inverse  of  the 
expunged  circulant  matrix  (if  the  inverse  exists)  does  not 
exist.  Using  the  optimum  detector  of  eq.  5.94  as  a 
guideline  we  hypothesize  that  the  following  detector  will 
have  near  optimum  performance.  We  base  our  decisions  on 

2(^)  where 
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N  ^  L-1 

:(X)  ^  y~ z(z|P,^)  -  y{n))* 

fer  ^  ^ 


N 

(1  +  ^a^Ll/p^_| 

1,. 


-  ^In 


N 

n  (1  +  A. 1  L/Pu  I ) 

k»1 


where 


P|,-)  Li./P|t_|-> 

•  4  w£ 


(5.98) 


-  '^k.lnd  +  A.l  L/P,.  1) 

K  —  H— 


-Co 


N 


VPk-l 


N 

'  Die'-’ 

k^T 


\  A|^|_^L/P|^^| 

m  1  *  AkLi*/^k-l 


VPw.l 


(5.99a) 


(5.99b) 


g  AijLt'/Pk-l 


(6.99c) 


For  each  PRS,  the  above  detector  utilizes  the  ciaximum  aumber 
of  integer  periods  available,  this  is  in  contrast  to  the 
optimum  detector  which  would  also  utilize  the  *'le£t  over” 
samples.  We  will  soon  show  that  when  L  is  large  the 
performance  of  the  above  sub-optimum  detector  is  nearly  the 
same  as  that  of  the  optimum  detector. 

m* 

it  ve  write  Lj^  for  Ll/P|j_|  the  form  for  z{^) 


reduces  to 


fig. 


5.  ■%  Averaged  R.O.C.  curves  for  the  quasi  opt xmu» 
detector  for  N-2,  each  PRS  has  identical  povef 
levels  of  .001,  .002,  .004,  .008  and  .016. 


[♦] 


I 
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.80  .91 


Averaged  R.O.C.  curves  for  the  quasi  sum 
detector  for  K">2.  with  the  signal  and  noise 
statistics  used  in  fig.  5.15. 


Pig.  5.17  R.O.C.  curves  for  the  energy  detector  with 
the  same  signal  and  noise  statistics  as 
used  in  figs.  5.15  t  5.16. 
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that  obtained  in  eg.  5.94.  We  call  the  detector  that  bases 

its  decisions  on  z{^)  the  quasi  optimuun  detector.  R.O.C. 
curves  for  the  quasi  optimum  detector,  for  N=2,  are  given  in 
fig.  5.15.  We  let  L*1000,  Pi  is  fixed  at  21  P2  is  allowed 
to  vary  from  22  to  30  and  the  R.O.C.  is  the  average  R.O.C. 

N 

We  call  the  detector  that  bases  its  decision  on  >  z^zlPjj) 

the  quasi  sum  detector.  The  R.O.C.  curves  for  the  quasi 
sum  detector  are  given  in  fig.  5.16.  It  is  seen  that  the 
performance  of  the  quasi  sum  detector  is  essentially 
equivalent  to  that  of  the  quasi  optimum  detector.  For 
comparison  purposes  the  R.O.Cw  curves  obtained  using  an 
energy  detector  are  given  in  fig.  5.17. 

An  approximate  expression  for  the  normal  detectability, 

d,  for  the  qua  si -'optimum  detector  is  given  in  eq.  5.100.  It 
is  assumed  that  N*2  and  that  both  the  PRS's  have  the  same 
power  level  A/2  and  that  Pi “Pa* 


d  a  A  L(1  +  Li  La) 


(5.100) 


The  normal  detectability,  ^or  the  optimum  detector  is 

then  approximately  upper  bounded  by 


<  A  L(1  +  Li  1  +  La  +  1) 
opt  ‘ 


(5.101) 


So  we  have 


I  “  (1  Li 
°opt 


L2)/(3  +  Li  ♦  La) 


(5.102a) 


in  the  limit  L  and  hence  Li  and  La  become  large 


1 


(5.102b) 


d 

V 

2 

For  the  energy  detector  we  had  (eq.  5.52)  cig^=A  L/2.  We 
compare  the  input  SNRs  for  fixed  performance  and  obtain 

I  =  (1  +  Li  +L2)^/^  (5.103a) 

Comparing  the  output  SNR  for  fixed  input  SNR  we  obtain 

I  »  1  +  Li  +  La  (5.103b) 

°ed 

The  above  results  are  analogous  to  the  results  obtained  in 
eq  5.54  and  eq.  5.55  and  are  in  good  agreement  with  the 
simulation  results. 

5.3.2  Detection  of  k  of  H  PRS* a  in  Noise 

This  is  a  generalization  of  the  detection  of  one  of  M 
PRS*s  problem  considered  in  Sect.  5.2.5.  Now  we  consider 
the  signal  to  be  the  sum  of  k  independent  PRS's«  where  the  k 
PRS's  belong  to  a  set  of  M  independent  PRS's.  We  assume 
that  it  is  not  known  a  priori  which  k  of  the  M  PRS's  make  up 

the  signal.  This  implies  there  are  N  -  (^)  possible  sets  of 
signals.  Even  for  moderate  M  and  k  the  signal  set  becomes 
quite  large.  As  we  will  show  the  optimum  detector  for  this 
situation  has  a  fairly  elaborate  and  complex  structure. 
This  fact  will  provide  us  with  the  motivation  for  searching 
for  detectors  which  have  performance  close  to  that  of  the 
optimum  detector  and  yet  have  much  simpler  structure. 

As  usual  we  will  assume  there  are  L  samples  in  the 
observation,  where  L  is  large  and  arbitrary.  The 
observation  under  the  hypotheses  Ho  and  is 
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Ho  :  y(i)  »  n(i)  ~  N(0,1)  i.i.d.  (5.104a) 

Hi  :  y(i)  =  n(i)  +  8^.(1)  »  s  &  n  independent  (5.105b) 

where  c  S  =  (si  ,82 , . . .  ,8^^} .  Each  of  the  elemente  of  S  is 

the  sum  of  k  PRS's  belonging  to  C  *  {pi  rfiar  •  •  •  •  The 

optimum  detection  statistic,  as  in  Sect.  5.2.5  is  the 
average  likelihood  ratio.  Since  L  is  not  necessarily 
composite  we  will  instead  use  the  average  of  the  conditional 
quasi-likelihood  ratios  as  the  "optimum"  detecton  statistic. 
Assuming  that  each  of  the  elements  of  S  is  equally  likely  to 
be  the  signal  we  obtain 


N 

S  y~exp(z(x|s^.) 


r»1 


(5.105a) 


N 

z(x)  ■  ln(^  exp(z(y|s^))  -  InN  (5.105b) 

To  develop  the  structure  of  the  quasi-optimum  detector  we 

need  to  expand  z(xtSj.).  Assume  s^  consists  of  PRS’s  of 
periods  »Pf2' •  *  *  »^rk'  ^^®  period  Pj.j^  has 

power  Apj^.  Now  we  can  write  z(x|Sr)  the  style  of 

eq.  5.97 


1  •jH 
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2(l|Sp)  •  z(z|Br^,Pr2»-*»Prk) 


L-1 

n«0 


-iln  ^ 

T  ^ - 

n  (1  +  Ap^Lpj^) 


1«1 


(5.106) 


where 


I-tl  ^  l-t'/Prl- 


rl- 


(5.107) 


L^x' 


'rl 


a  A  A-i 

or  - 1 


n 

'•£ 


I 


^rX^rl 


^rl^rl 

•T  ^  ^rl^rl 


(5.108) 


P,l-1  L,i-1 


2<Z|Pri)  •  <y” y<j+sPpl)) 


^  3=^  8=^ 


-  ^rl.lnd  + 


(5.109) 


and 


c  1  A  ^rl^rl 

TTir-X 


rl“rl 


(5.110) 


The  detector  first  calculates  the  M  conditional  quasi  log- 


likelihood  ratios  zC^lPi)#...,  Then  it  forms  (|^) 

sums  of  z(^|Pj)  in  groups  of  k  and  finally  adds  the 
previously  calculated  terms  to  the  sums  to  obtain  z(y|Sj.), 
The  terms  z(ylsi z(y|Sjj)  are  exponentiated  end 

m0 

averaged  to  obtain 

The  estimator-detector  bases  its  decisions  on 
max  z(y|s  ).  The  structure  of  this  detector  is  not  much 

Sj,cS 

M 

simpler.  In  fact  it  also  caculates  N  =  (j^)  terms  and  bases 
its  decision  on  them. 

In  Sect.  5.3.1  we  had  seen  that  basing  decisions  on 
N_ 

z(^|Pj^)  was  essentially  equivalent  to  basing  decisions  on 


z(y).  The  next  detector  we  try  bases  its  decisions  on 
k 

max  r;  (vlP  -.).  But  this  is  the  same  as  basino  the 
r«1 , . . .  ,N  “•  tA 

decision  on  the  sum  of  the  k  largest  z(y;|P^),  i«1,...,M. 

This  detector  has  a  much  simpler  structure  than  the  previous 
two.  We  call  this  detector  the  max-detector.  Choosing  the 

k  largest  z(y|P>}  also  immediately  forms  an  estimate  of  the 
periods  of  the  k  PRS's  in  the  signal. 

The  R.O.C.  curves  for  the  quasi-optimum  detector,  the 
estimator-detector  and  the  max-detector  are  given  in 
figs.  5.18,  5.19  &  5.20  respectively.  We  let  M-10,  k»2  and 
n«{21 ,22, . . . ,30} .  We  see  that  the  performance  of  the  max- 
detector  is  nearly  identical  to  the  performance  of  the 
estimator-detector  but  slightly  poorer  than  the  performance 
of  the  quasi -optimum  detector. 


Fig  5.19  Averaged  R.O.C.  curves  for  the  estimator 
detector  for  M>10  and  k»2  with  the  same 
signal  and  noise  statistics  as  in  fig.  5.18. 
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I 


Fig  5.20  Averaged  K.O.C.  curves  for  the  max  detector 
for  M"10  and  k«2.with  the  same  signal  and 
noise  statistics  as  in  figs.  5.18  I  5.19. 
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5.3.3  Detection  of  Arbitrary  PRS  in  non  White  Gaussian  Noise 

To  complete  the  discussion  on  the  detection  of  PRS's, 
we  examine  the  situation  where  the  signal  is  complex  and 
successive  signal  samples  are  not  independent,  we  also 
assume  that  the  additive  noise  is  not  necessarily  white.  We 
assume  that  the  signal  is  at  least  wide  sense  stationary  and 
that  its  autocorrelation  matrix  is  circulhnt.  Signal  and 
noise  are  still  assumed  to  be  mutually  independent.  We  will 
derive  the  equations  and  give  conditions  under  which  they 
are  valid  without  going  in  great  detail.  We  assume  there 
are  L  samples  in  the  observation,  where  L  is  some  integer 
multiple  of  the  signal  period.  Let  R^^  be  the  noise 
autocorrelation  matrix  and  R^  be  the  signal  autocorrelation 
matrix.  The  observation  statistics  under  hypotheses  Hq  and 
Hi  are 

-  (j;)'''^^|Rnr’''^exp(- ^2*r;;’2)  (5.111a) 

£(2|Ki)  ■  (5.111b) 

From  the  above  we  may  write  down  the  likelihood  ratio  and 
the  log-’likelihood  ratio 

-  iRyjl |Rj^+Rgr'/^.exp(^y*(R’^  -(Rj^+Rg)"^^)  (5. 1 12a) 
z(^)  -  “  3^nl^n'"^sl  (5.112b) 

In  the  above  equations  we  have  assumed  that  the  appropriate 
inverses  exist.  Rg  is  circulant  due  to  the  signal  being 
periodic.  If  the  noise  is  white  R^  will  be  the  identity 
matrix  and  hence  trivially  circulant  and  also  invertible. 
If  the  noise  is  not  white  R^  will  be  circulant  if  the 
associated  autocorrelation  sequence  R^Ot)  satisfies  the 
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following  condition 


Rn(n)  »  R*(L-n)  (5.113a) 

or 

Rn(n)  «  Rj^(L+n)  (5.113b) 

In  other  words,  if  the  noise  is  not  white  we  require  that 

the  noise  and  the  noise  autocorrelation  be  periodic,  and  the 

number  of  samples  in  the  observation  L  either  be  the  noise 

period  or  some  integer  multiple  of  the  noise  period.  With 

R„  and  R-  circulant  we  have 
n  s 


(5.114.) 


R  ■  u-'f  ur 


(5.114b) 


R.+R-  •  U,  (y.+i(.)o; 

tl  O  M  II  0 


(5,114c) 

(5.114d) 


(R  +R  )*^  ■  Ur(f  )“^c5 

'  n  s'  L'^n  ^s' 


(S.114e) 


where  and  Yg  are  diagonal  matrices  and  the  diagonal 
elements  are  given  by  t(ig(m)  and  ^g(m)  respectively 


L-1 

♦„(m)  . 


(k).exp(-  j2wkm/L) 


(5.115a) 


L-1 


-  jy, 

Ic^ 


4»«(ro) 


(k).exp(-  j2irkm/L) 


(5.115b) 
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Now  for  r”^  to  exist  we  require  that  >  0  for  all 

m,  this  means  that  Rfj(h)  can  not  have  a  period  smaller  than 

L.  Whenever  >  0  for  all  m  both  rJJ^  and  yiill 

exist.  Substituting  for  r"^  and  ®Q*  5.112b  we 

obtain 


z(^) 


I'^n' 


('F  +'{'-) ’bu*y  -  ilnl'^n'^’^s 
ns 


(5.116) 


now 


Z  \  -  l"’/^.(7(0),Y(-1),...,?(-L+1))  (5.117a) 

where  Y  is  the  complex  conjugate  of  Y  also 

U*x  -  L“^/^.(Y(0),Y(-1),...,Y(-L+1))'^  (5.117b) 

Now  let  -  (Yj^-^-Yg)"^  ^  0,  0  is  a  diagonal  matrix  and  let 
the  diagonal  entries  be  6(m)  then  z(^)  is  given  by 


L-1 


y*~0(ro)  |Y(~ro)  I  ^  - 


m«0 


L-1 

^  e(m)|Y(L-m)l* 

m*>0 


ainlV^s! 

^  ~IV 


(5.118a) 


(5.118b) 


The  above  detector  may  easily  be  generalized  for  all  the 
cases  discussed  in  the  earlier  sections  using  procedures 
similar  to  the  ones  used  earlier.  The  above  analysis  holds 

for  white  noise  if  we  replace  R^,  r”\  Y^^  and  y”^  by  the  LxL 
identity  matrix. 
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5»4  The  Abstract  Problem  Revisited 

In  Sect.  5.1  we  had  developed  the  analogy  between  the 
exact  problem  developed  in  Chapter  III  and  the  detection  of 
two  independent  PRS's  pi(n)  and  psCn)  with  relatively  prime 
periods  P^  and  Pa  respectively  and  an  observation  interval 
of  L-P^.Pa  samples.  Using  this  analogy  and  the  results 
established  in  Sects,  5.2  &  5.3  we  will  answer  the  questions 
asked  at  the  begining  of  this  chapter  and  establish  some 
guidelines  for  the  exact  problem. 

For  the  two  PRS  problem  we  considered  two  detection 
philosophies  1)  Consider  the  L^Pj.Pa  samples  to  be 
independent.  2)  Consider  the  L  samples  as  having  been 
generated  from  Pi^^Pa  independent  samples.  According  to  the 
first  philosophy  the  optimum  detector  was  the  energy 
detector,  the  second  philosophy  led  to  the  circular 
averaging  energy  detector  (frequency  selective  energy 
detector)  as  the  optimum  detector.  The  performance  of  the 
CAEO/FSBD  was  Shown  to  be  considerably  better  than  that  of 
the  (ci9S.  o.ii  a  S.'i3  and  eqs.  5.102  4  5.105).  This 
suggests  that  the  detector  for  the  original  problem  have  N'*>H 
degrees  of  freedom;  that  is,  it  should  consider  the  HH 
peaks  in  the  ambiguity  diagram  as  having  been  generated  from 
N-^N  independent  sets  of  doppler  and  delay  parameters.  As  an 
example  let  H«9  and  H**?,  i.e.  we  have  9  paths  to  one 
receiver  and  7  paths  to  the  other  receiver.  The 
crosscorrelation  of  the  receptions  at  the  two  receivers  will 
have  63  peaks.  Treating  the  63  peaks  as  being  generated 
from  16  independent  parameter  sets  rather  than  as  63 
independent  entities  results  in  a  performance  improvement  of 
the  order  of  12dB  (obtained  by  substituting  for  and  L2  in 
eq.  5.103b). 

The  near  optimum  performance  of  the  estimator "detector 
and  the  max"detector  (figs.  5.16,  5.19  t  5.20)  suggests  that 
the  estimates  of  the  path  parameters  (differential  doppler 
and  delay)  be  used  to  locate  the  the  MN  peaks  in  the 


crosscorrelation  or  the  ambiguity  diagram.  Differential 
doppler  and  delay  estimation  is  a  topic  of  on-going  research 
interest  and  is  outside  the  scope  of  this  thesis. 
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CHAPTER  VI 

SUMMARY  AND  CONCLUSIONS 

The  object  of  this  thesis  was  to  form  guidelines  and 
rules  of  thumb  for  an  experimenter  faced  with  the  problem  of 
detecting  a  moving  acoustic  source,  where  sound  is  assumed 
to  propagate  over  a  multipath  channel.  Source  signals  were 
received  at  two  remote  receivers;  the  Bivariate  Normalized 
Crosscorrelation  (BNC)  function  (and  the  associated 
ambiguity  function)  of  the  receptions  at  the  two  receivers 
was  used  as  the  observation  for  the  detector.  It  was 
assumed  that  there  were  M  paths  to  one  receiver  and  N  paths 
to  the  other  receiver  and  that  the  observation  times  were 
long  enough  for  the  path  pairs  to  be  resolved. 

The  exact  multipath  channel  problem  is  difficult  to 
solve,  as  was  shown  in  Chapter  111,  so  most  researchers 
asssume  that  the  propagation  is  by  a  single  path  channel. 
Detectors  based  on  the  single  path  channel  model,  when  the 
propagation  is  by  multipath,  are  by  definition  sub'‘optimum 
because  they  do  not  utilize  oil  the  information  available 
about  the  source  in  the  receptions.  Our  approach  was  novel 

in  the  sense  that  we  chose  to  retain  the  multipath  channel 

assumption.  Based  on  this  we  formulated  two  solvable 
abstract  problems  for  further  study.  The  abstract  problems 

were  chosen  to  be  sufficiently  close  to  reality  so  that 

their  study  would  provide  useful  information  about  the  real 
problem.  This  also  provided  a  framework  for  maximizing  the 
information  available  to  the  detector  about  the  source. 

In  Chapter  IV  we  studied  the  "Extended  M-‘Orthogonal 
Signals"  problem.  This  provided  us  a  method  for  combining 
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the  information,  available  in  the  ambiguity  plane,  for 
detection  purposes.  We  also  showed  that  it  was  necessary  to 
localize  a  quarter  to  a  half  of  the  peaks  in  the  ambiguity 
plane  for  acceptable  detection  performance. 

The  ambiguity  diagram  consists  of  MN  peaks  due  to  the 
signal.  The  MN  peaks  are  generated  by  independent  sets 
of  doppler  and  delay  parameters.  In  Chapter  V  we  addressed 
the  question;  Should  the  detector  treat  the  MN  peaks  as 
independent  entities  or  as  being  generated  from  M-i-N 
independent  parameter  sets?  In  other  words,  should  the 
detector  have  MN  degrees  of  freedom  or  M+N  degrees  of 
freedom?  To  answer  this  question  we  used  Periodic  Random 
Sequences  (PRS's).  It  was  shown  that  a  detector  with  M-^N 
degrees  of  freedom  had  better  detection  performance  than  a 
detector  with  MN  degrees  of  freedom.  Also  the  results  of 
Chapter  V  suggest  that  the  M-»-H  sets  of  doppler  and  delay 
parameters  be  estimated  and  the  estimates  be  used  to 
localize  the  peaks  in  the  ambiguity  diagram. 

Putting  the  results  of  Chapter  IV  and  Chapter  V 
together  we  can  form  the  following  guidelines: 

1)  Use  observation  (or  integration)  times  that  are  long 
enough  so  that  when  the  BNC  function  is  formed  path  pairs 
can  be  resolved. 

2)  Estimate  the  M'^’N  sets  of  doppler  and  delay 
paramters. 

3)  Use  the  estimates  to  form  the  differential  doppler 
and  differential  delay  values,  and  use  these  values  to 
localize  the  peaks  on  the  ambiguity  diagram, 

4)  Use  the  Extended  M*Orthogonal  Signals  technique  to 
combine  the  peak  values  in  the  ambiguity  diagram. 

Acceptable  detection  performance  is  obtained  if  at 
least  a  quarter  to  a  half  of  the  peaks  have  been  correctly 
localized. 

The  Extended  M-Orthogonal  signals  formulation  is  based 
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on  the  classical  M-Orthogonal  signals  problem,  however  it  is 
a  new  formulation.  The  theory  developed  for  the  detection 
and  estimation  of  Periodic  Random  Sequences  is  based  on  well 
known  results  of  detection  and  estimation  theory  but  this 
formulation  is  believed  to  be  new.  This  has  application 
wherever  the  signals  exhibit  periodic  or  cyclic  behavior 
(engine  noise,  cyclic  codes  etc.). 

Further  work  is  needed  in  passive  estimation  of  doppler 
and  delay  parameters  assuming  multipath  propagation.  This 
could  form  the  basis  for  a  Doctoral  Thesis.  Given  the 
estimates  of  the  doppler  and  delay  parameters,  localization 
of  the  peaks  in  the  ambiguity  plane  could  form  the  basis  for 
a  masters  project.  To  test  the  validity  of  the  guidelines 
established  in  the  thesis  (and  their  refinement) 
experimental  work  needs  to  be  done. 
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APPENDIX  A 

MONTE  CARLO  SIMULATION  METHODS 

The  R.O.C.  curves  and  the  scatter  plots  in  Chapter  V 
were  generated  using  Monte  Carlo  simulations.  In  this 
Appendix  we  use  the  example  of  one  PRS  in  noise  to  briefly 
outline  the  simulation  methods. 

Assume  that  the  observation  interval  is  L  samples  long. 
Let  the  PRS  period  be  Pi  and  the  PRS  power  level  be  Ai.  For 
convenience  we  let  L»LiPi,  where  Li  is  the  number  of  integer 
periods  of  the  PRS  in  the  observation.  We  repeat  the 
equation  for  z(y),  the  detection  statistic,  in  modified  form 
below 


PrlLi>l 

2(Z)  »  Ki  y~ (Y" y(j+kPi)>*  -  Kj  (AI) 

3^  k'«d 

Ki  and  Kj  are  constants  that  do  not  depend  on  the  reception 
and  can  be  pre-calculated. 

The  probabilities  of  detection  and  false  alarm  were 
calculated  on  the  basis  of  2,000  independent  trials  under 
both  Hi  and  Ho*  Let  us  assume  that  the  vectors  vo  and  Vi 
each  contain  2,000  independent  outcomes  of  z(y)  under  Ho  and 
Hi  respectively.  Two  25  bin  histograms  HGo  and  HGi  were 
formed  using  the  data  from  Vq  yi  respectively.  The  25 
bins  for  both  the  histograms  were  "uniformly  spaced"  between 
min  yi  and  max  vo*  Data  points  from  vq  smaller  than  min  yi 
were  put  in  the  first  bin  of  HGo.  Data  points  from  Vi 
larger  than  max  yo  were  put  in  the  last  bin  of  HGi.  With 
the  25  bin  histograms  available  it  is  fairly  straightforward 
to  calculate  25  points  on  the  R.O.C.  curve. 

Now  we  describe  the  generation  of  the  vectors  Vo  and 
Vi.  A  uniform  random  number  generator  that  generated 
independent  variates  in  the  range  [0,1)  (U  "  (0,1))  with  a 
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.  2* 

period  of  2  -1  was  available.  The  Box-Muller  method.  [38] 

was  used  to  generate  zero  mean  unit  variance  gaussian 
(N~(0,1))  random  variates.  Given  M  i.i.d.  U  ~  (0,1)  random 
variates  the  Box-Muller  method  gives  a  simple  and  efficient 
procedure  for  generating  M  i.i.d.  N  "  (0,1)  random  variates 
(where  M  is  even).  The  elements  of  Vo  were  generated  as 
follows 

Step  1  :  Generate  L  uniform  random  var5ates  if  L  is 
even  L+1  if  L  is  odd. 

Step  2  :  Use  the  Box-Muller  method  io  get  L  N."  (0,1) 
random  variates. 

Step  3  :  Perform  the  signal  processing  indicated  in 
eg.  A1  to  get  an  element  of  Vo. 

Step  4  :  Repeat  steps  1  through  3  until  Vq  contains 
2,000  elements. 

The  elements  of  Vj  were  generated  as  follows: 

Step  1  :  Generate  U  uniform  random  variates  if  L  is 
even,  L+i  if  L  is  odd. 

Step  2  :  Use  the  Box-Mullet  method  to  get  L  N  *  (0,1) 
random  variates  and  store  in  a  vector. 

Step  3  :  Generate  Pi  uniform  random  variates  if  Pi  is 
even,  Pi+1  if  Pi  is  odd. 

Step  4  :  Use  the  Box-Muller  method  to  get  Pi  N  "  (0,1) 
random  variates. 

Step  5  :  Multiply  the  Pi  N  **  (0,1)  random  variates  by 

step  6  :  Periodically  extend  the  Pi  random  variates  in 
Step  5  Li  times  and  store  in  a  vector. 

Step  7  :  Add  the  vectors  of  Step  2  and  Step  6  element 
by  element. 

Step  8  :  Perform  the  signal  processing  indicated  in 
eg.  A1  to  get  an  element  of  vi. 


step  9  :  Repeat  steps  1  through  8  until  Vj  contains 
2,000  elements. 

The  above  procedures  give  an  idea  of  the  simulation 
methods  used  for  the  more  complicated  problems. 
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DERIVATION  OF  THE  NORMAL  DETECTABILITY  OF  THE  FSED 


In  this  appendix  an  expression  for  the  normal 
detectability  of  the  FSED  (CAED)  is  derived.  We  assume  that 
the  signal  consisists  of  two  PRS’s  with  periods  Pi  and  P2. 
We  let  the  observation  be  L  samples  long,  where  L=Pi.P2.  We 
define  L|^=L/P|^  to  be  the  number  of  periods  in  the 
obsesrvation  of  the  PRS  with  period  Pj^.  We  further  assume 
that  Pi=iP2«  The  normal  detectability  is  defined  as 


[Ei(2)  -  Eo{z)]' 
var{z) 


(B1) 


We  will  use  the  variance  of  z  under  Ho  in  eg.  B1. 
Expresssions  for  z(^)  for  the  situation  described  above  are 
given  in  egs.  5.78b,  5.86,  5.94  and  5.95.  Using  egs.  5.87, 
5.88,  5.89,  5.90  and  5.91  an  alternative  expression  for  z(;]^) 
is  obtained 


z{^)  «  {2L) 


L-1 

% 

m»5 


0(m) |y(m) 


Inl'Fl 


{B2) 


The  last  term  in  eg,  B2  does  not  affect  the  value  of  d  and 
will  be  carried  through  as  a  constant,  ’’K”.  Y(m)  is  the  L 
point  DFT  of  the  observation 


Ho  :  Y(m)  -  N(m) 


(B3a) 


Hi  :  Y(m)  «  N(m)  *  Pi(m)  h  p2(m} 


(B3b} 


where  N(m),  Pi(m)  and  P2(in)  are  the  DFT's  of  the  noise,  PRS 
with  period  Pi  and  PRS  with  period  P2  respectively.  6(m)  is 
defined  as  (see  egs.  5. 90  and  5.91} 
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Pi-1 


P2-I 


e(m)  =  “ 


Co6(m)  +  Cl  ^  6(m“rLi)  +  C2  /  6(m-rL2) 


(B4) 


Eo(z)  =  (2L) 


L-1 

0(m)Eo(|y(m)  1^ 


(BS) 


Expanding  Y(m)  we  obtain 


L-1  L-1 


Eo(|Y(m)|  ) 


Bn(k)n(l)  .exp(  j2iTm(l-k)/L)  (B6a) 


L-1  ^ 


6(k-l)  .exp(  j2‘<rm(l-k)/L)  (B6b) 


it  follows  that 


Eo(z)  •  j  (“Co  CiPi  C2P2)  -  K 


(B7) 


Bo(z)  -  (2L)”^  y”  L.e(m)  -  K 


(B8) 


0(m)  ■  -Co  Cl?!  +  C2P2 


(B9) 


(BIO) 


Now  we  calculate  the  expectation  of  z{^)  under  Hj. 


Ex(2) 


L-1 

(2L)’^  ^ 


0(m)Ei(|Y(m)|  )  -  K 


m 


(B11) 


the  expectation  within  the  summation  can  be  written  as 
Ei(lY(m)|^  =  E(|N(m)  +  Pi(m)  +P2(m)|*) 

-  E(|N(m)|^  +  |Pi(m)i*  +  |P2(m)|*  (812) 

+  2Re(N(m)pt(m)  +  N(m)pJ(m)  +  Pi(m)P?(m))) 

Since  the  PRS’s  and  noise  are  mutually  independent  and  zero 
mean  the  above  reduces  to 

Ei(lY(m)l^  -  E{|N(m)|*)  +  E(|Pi|*)  +  E{|P2(m)|*)  (B13) 

2 

We  already  know  E(|N(m)|  ),  we  only  need  to  find 
2  2 

E(|Pi(m)|  ),  E(lP2(m)|  )  will  have  a  similar  expression. 

Now 


L-1 

Pi<m)  ■  Pi  (1)  ,exp(-j2irlm/L) 


Li-IPi-1 


Pi(s)  .exp(-j2wsffi/L)  .exp(-j27rkm/Li)  (B14) 


where  pi(s)  is  one  period  of  pi(l).  Eq.  BU  simplifies  to 
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Pi-1 


Pi{!n)  «  Li  )  Pi  (s)  .exp(-j2irsv/Pi ) 


(B15) 


m  »  vLi  ,  V  ■  0,1,.. . ,Pi-1 


It  follows  that 


Pi-IPi-1 


|Pi(m)|  »  Li  2  2  Pi(s)pi(k).exp(j2irv(k-s)/Pi) 

s«0  k»0 

m  »  vLi  ,  V  ■  0,1,...,Pi-1  (B16) 


and 


Pi-IPi-1 


B  Pi  (in)  I  ■ 


Epi(s)pi  (k)  .exp(  j2irv(k-8)/Pi) 


Pi-1Px-1 


2 

AiLi 


6(k-8)  .exp(  j2Trv{k“S)/Pi>  (B17) 


8« 


m  ■  vLi  ,  V  ■  0,1,...,Pi-1 


The  above  simplifiea  to 


Pi~1 


E(|Pi(m)|  )  ■  AiLiPi  )  6(o-vLi) 


JPi  Y" 

v»0 


(B18) 


similarly  we  have 


Pa-1 


E(|Pi(m)|*)  -  AaL^Pa  Y~  ^(m-vLa) 

v*0 


(B19) 


V\''V  .  .  ; 


.•  .*•  ..'V  *.•  W.*  >  V  V.’.’ 

■  ■  .•  <L*  «>■  A  .•>  .*>  .v.*>  V f.'- 


2  T — 

-  AiLjPi  ) 

v«0 


(-Co6(vLi)  +  Cl 


6( (v-r)Li) 


P2-I 

+  C2  ^  6(vLi-rLs)) 

r»0 

2 

•  AiLiPi(—Co  ^  *  C2) 


Similarly 


L-1 

3)  AaWPa 

■  AjL^t-Cft  ♦  Cl  ♦  C2P2) 

It  follows  that 

£1(2)  •  (“Co+CiPi+CjPa)  (2L)  .AiLiPiC-Co+CiPi-^Cj) 

♦  (2L)“^A2L5Pa(-Co’K:i+CaPa)  -  K 

■  4  (-Co  CiPi  ♦  CaPa  ♦  AiLi(-‘Co  CiP|  Ca) 

♦  AaLa(-Co  ♦  Ci  ^  CaPa))  “ 


XT  XT 


m»0  v*0 


We  can  now  write  an  expression  for  the  difference  in 
expectations  of  z  under  Hi  and  Ho 


(B24) 


(B2S) 

the 


Ei(z)  -  EoCz)  “  ^  (AiLi(C2-Co)  +  A2L2(Ci-Co) 

+  CjAiL  +  C2A2I1) 


Now  W8  find  the  variance  of  z  under  Ho 

2  2 

varo(z)  *  Eo(z)  -  (Eoz) 

2 

For  z  we  have  the  following  expression 


2 

z 


L-1  L-1 


6(m)e(n)|y(m)|^iY(n)|^ 


L-1 

y~  e(m)  ly(m)  +  K* 
maO 

2 

Eo(z)  is  then  given  by 

L-1  L-1 

Eo(z)^  a  (2L)"^  Y~  9(m)e(n)Eo|Y(m)  |*|y(n)  I 

fej' 

L-1 

“I  y~  9(*n)EolY(m)  I*  +  K* 
m»(y 

2 

and  (Eoz)  is  given  by 
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